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ABSTRACT 


Training residents for orthopedic surgery can be challenging due to the 
variety of parameters involved in critical tasks. This project presents a method 
for skill assessment of surgical drilling. The parameters measured, which include 
arm and tool kinematics, applied force, tool and bone vibration, and drill RPM, 
were measured using a combination of force, acceleration, and optical tracking 
sensors. High correlations were found between different combinations of pene¬ 
tration distance, applied force, drill toggle, drill orientation, and drill RPM. As 
a whole, the differences between experienced and novice surgeons measured in 
these trials were not statistically significant. However, when looking at certain 
performance criteria individually, "expert" surgeons tended to outperform their 
"novice" counterparts. As the second part of the study, we explored the feasibility 
of using a predictive model of axial force during simulated orthopaedic drilling. 
Synthetic bone samples were used with a robotic device to simultaneously record 
force, position, and velocity trajectories. A nonlinear function relating force to ax¬ 
ial position and velocity was used to fit the data. The models obtained for various 
bone types fit the well, and are relatively distinct, one from another. This suggests 
that a predictive model can be used and is able to capture relevant variations in 
the thickness and hardness of cortical and cancellous bone. The device used to 
capture this data may also be utilized for augmented haptics, allowing the sim¬ 
ulation of anatomic variability in bone while providing necessary force/torque 
levels. 

Keywords: Bone drilling; Skill assessment; Augmented haptics; Orthopaedics; 
Virtual reality 
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CHAPTER 1 


INTRODUCTION AND BACKGROUND 


1.1 Motivation 

Standard tools from the field of robotics may allow significant advances 
in surgical-skill assessment, training, and related procedures. In 2009, the Pa¬ 
tient Safety Committee of the American Academy of Orthopaedic Surgeons re¬ 
leased results of a member survey, which was used to identify common errors in 
orthopaedic surgery [3]. Results of this study showed that the orthopaedic sur¬ 
geons were directly involved in 60% of the total errors. Seventy eight percent of 
errors occurred in the hospital, with 54% occurring in the operating room. The 
hand and/or fingers and the knee (35% for each) were the most frequent sites of 
surgical error. As of July 2013, prompted by the quality initiative set forth by the 
Patient Protection and Affordable Care Act, the American Board of Orthopaedic 
Surgeons (ABOS) and the Residency Review Committee (RRC) for Orthopaedic 
Surgery have mandated formal motor skills training outside of the OR [1, 6] Fig¬ 
ure 1.1. The ABOS and RRC believe that this training will enhance resident's 
acquisition of basic surgical skills. Currently, no criteria has been established 
for hands-on training of surgical residents on common procedures used in each 
orthopaedic surgery subspecialty rotation. Since the inception of the ABOS re¬ 
quirement, some surgical skills simulation has been implemented by residency 
programs in the form of fracture skills courses [6, 2, 7], Unfortunately, limited 
data exist on the effectiveness of these courses. Prior work in surgical skill eval¬ 
uation using sophisticated robotic training tools includes the work specifically 
where data obtained through the actual surgical robot [8, 9] or by video during 
the surgery [5] is subsequently analyzed. This study used various robotic tools to 
improve the quality of training in different fields. 


1.2 Background and previous work 

1.2.1 Quantifying and Improving Surgical Results in Orthopaedics 

Although modern teaching skills mentioned above are important aspects 
of surgical training, quantified assessment has not played a significant role in it. 
Orthopaedic surgical residents were traditionally taught in operating rooms were 
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they observe surgical tasks. Currently, skill laboratories are also being used to 
teach and test residents. These skill laboratories not only improve the confidence 
of the residents significantly but also improve the results in orthopaedic surgery 
[10]. In general, research in this field follows two main goals: improving the 
tools in the operation room (OR) and improving training methods by quantifying 
results. 



Figure 1.1: A robotic system to detect and prevent over-penetration is shown[l]. 


1.2.1.1 Improving the tools in the operation room 

In last three decades, a wide variety of approaches have been applied to 
help surgeons in the OR ranging from developing simple tools like drill guides 
to sophisticated, robotic tools. Although autonomous drilling robots might have 
high precision, using these robots in ORs is not yet possible due to safety fac¬ 
tors [1]. These robots can monitor many critical parameters during the task, but 
variations in shapes and properties of human organs also limits the automations 
capacities [11]. 

Numerous skills are critical to orthopaedic repair. Over-penetration is one 
of the biggest challenges in orthopaedic drilling. To detect the over-penetrarion, 
multiple studies with different algorithms were conducted [1, 12], and different 
automated robots used some of the results in the mentioned studies to predict 
this action. Some of these studies were general and some others were for specific 
operations, such as spine surgery. Figure 1.2 [2]. Although most of these robots 
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monitor the revolution per minute (RPM) or the axial force, in some cases they 
apply image processing to find the location of the drill bit [6,13,14]. 



Figure 1.2: Optical methods to register the robotic arm for bone drilling[2]. 


Temperature is another critical parameter in bone drilling receiving signif¬ 
icant attention from scientists [15]. Overheating can cause serious, and in some 
cases permanent, damage to the bone and tissues structure in the vicinity Distinct 
experimental and theoretical approaches were applied to find the temperature of 
the bone during drilling. In general, experimental studies related to bone drilling 
show better results than numerical studies [16,17,18,19, 20]. Automated robots 
used most of these results and methods to perform a better task. 

Optimization of drilling has a significant role in bone drilling studies [19]. 
Unlike temperature and over-penetration, numerical and theoretical approaches 
to optimization have not yielded consistent results [7], Neural Network and 
Fuzzy logic have shown the best results compared to others [9]. Drill bit geom¬ 
etry, bone type, RPM, temperature, and down ward force are the most impor¬ 
tant parameters of the optimization^, 8, 9, 21, 19]. Combining the results of the 
optimization studies with safety studies is another challenge that the automated 
robots have to consider. The most optimized performance is not always the safest 
one. In general, robotic tools in the OR must follow various algorithms to be able 
to keep the task safe and optimized. 


1.2.1.2 Improving the training methods by quantifying the results 

Perhaps the most obvious approach to improving bone drilling in the op¬ 
erating room is to improve training methods for orthopaedic residents. One of 
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the main requirements in a precision training method is quantifying the test. To 
create the quantified scale and test, results from safety and optimization studies 
were combined with the skill assessment and monitoring approaches[22]. Com¬ 
paring experts with novice surgeons is an effective way to design a quantified 
scale. Figure 1.3[3,5]. For example, a quantified scale has been made by following 
the motion of the hand during the laparoscopic surgery and comparing certified 
and non-certified surgeons [4], Robb RA, developed a virtual reality system for 
endoscopy As a skill assessment, they developed specific clinical protocols that 
compared virtual endoscopy with real endoscopy [23]. 




(a) 


(b) 


Figure 1.3: (a): Using electro magnetic sensor to track the hand's motion[3]. (b): 
Using gyroscope to track the hand's motion[4]. 


There have been multiple studies to simulate the surgical tasks without 
preliminary output. Later, Lee et al. made another simulator for laproscopic 
surgery with a quantified skill assessment approach [24], In 2004, Eversbusch 
and Grantcharov found the learning curves for a simulated colonoscopy via a 
randomized trial algorithm [25]. Recently, Sezek et al. conducted research in 
training, simulation, and reducing complications in Urology. The study showed 
e-learning and observation of a procedure have a significant effect on the learn¬ 
ing curve in reducing the learning time and improved the learning peak[26]. Ue- 
mura et al. showed that, the number of operations performed by a surgeon do 
not necessarily improve the skill of the surgeon. Figure 1.4 [5]. Reznick applied 
educational psychology methods to optimize the of teaching residents [3]. 
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Figure 1.4: Comparing the movement of two hands for experts (left hand, left-up, 
right hand, right-up) vs. novice (left hand, left-down, right hand, right-down) [5]. 


Finally, tools can and are being designed with safety features to add ad¬ 
ditional safety measures to the procedure. In one case hand-held drills with an 
interruption mechanism were used to avoid over-penetration [27], This study 
showed the design was effective, but additional work is required to allay safety 
concerns. Designing a safe device can be a significant challenge given the uncer¬ 
tainty in human behavior. 


1.2.2 Development of Robotic Tools for Orthopaedics Surgical Skill Train¬ 
ing 

1.2.2.1 Haptics 

Haptic technology is designed to interact with the human touch sense. 
Quickly, this technology was incorporated into virtual reality (VR) and train¬ 
ing simulators. It also can play a significant role in teleoperation systems such 
as those used for explosive ordinance disposal or telesurgery [28, 29]. These 
robots can be classified into two major categories admittance and impedance 
types. Impedance type robots are most common in commercial haptic robots. 
Impedance robots are designed to be back-drivable by minimal friction and in¬ 
ertia. Admittance type robots can generate more force feedback but due to their 
high friction and inertia, they are not back-drivable [30]. 
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1.2.2.2 Augmented Haptics 


Although haptic robots are popular in VRs and training simulators, oc¬ 
casionally, they suffer from lack of fidelity with respect to graphical and haptic 
elements. The use of augmented haptics can largely address this problem. In 
2014, Steyn et al. developed an augmented-reality environment for a wheeled 
mobility system [31]. One of the main usage of haptic robots is in tele-operations. 
Maddahi et al., conducted a study to reduce the error of manipulator robots, by 
combining the haptic environment with a robotic assisted-line tracker[32], 

Augmented reality environments (AEs) do not always include haptics. For 
example, Hallet et al. combined the graphical VR with the real environment in the 
operation room (OR) to decrease the risk of error during the trans-thoracic (min¬ 
imally invasive) liver resection[33]. In contrast, there have been similar studies 
where haptics were included in the AE to great effect, such as haptic systems 
with an augmented observer approach for high-impedance haptic system, and 
laser projection with surface recognition for laparoscopic surgery [34, 35, 36]. 


1.2.2.3 Surgical Simulators 

Using robots to improve training methods is not a new idea. Simulators 
have been used in different studies to train surgeons. For simulators to be real¬ 
istic, a precision mechanical model is required[37, 19, 38, 39, 40]. Some of these 
models were derived from other studies in safety and optimization fields. Some 
of these models are flexible[41] while others are fixed, but with feedback for learn¬ 
ing purposes[42]. Even with a precision mechanical model, most of the haptic de¬ 
vices are not strong enough to simulate the exact feeling of the bone drilling[43]. 
Instead, the focus of the simulators are on graphics, sound, and feedback[44, 45, 
42]. The efficiencies of the feedback from these studies are sometimes unreliable, 
due to the lack of accurate force feedback in the simulators [11,43,41,42], Despite 
these limitations residents with simulator training were shown to preform better 
in surgical tasks and had higher self-confidence [46, 47, 48]. 


1.3 Organization of the Thesis 

The chapters that follow are organized into three major parts, which each 
showcase how tools from the field of robotics can be applied to current clinical 
needs: surgical-skill assessment (chapter 2), augmented haptics for bone drilling 
simulators (chapter 3). Finally, the fourth chapter contains the conclusion and the 
future recommendations. 
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CHAPTER 2 


ASSESSMENT OF ORTHOPAEDIC SURGICAL SKILL 


2.1 Introduction 

Orthopaedic training is a hands-on, motor skills-demanding surgical spe¬ 
cialty. Historically, much of the technical learning occurs in the operating room 
(OR). This is an effective method of training for simple fractures or low-risk cases. 
For more complicated operations or high-risk situations (spinal surgery, potential 
for vascular or nerve impingement, etc.), using the OR as an introduction to sur¬ 
gical techniques can introduce unnecessary risk. Instead the motor skills tasks 
can be simulated in the lab prior to performing the tasks in the OR. This model 
has been utilized in the endoscopic training of general surgeons. 

As was mentioned in chapter one, a standardized method of training out¬ 
side of the operating room has not been developed. Modules delineating the 
tasks deemed important by the ABOS have been provided as guidance in this 
endeavor and are publicly available (www.ABOS.org). Since the inception of the 
ABOS requirement, some surgical skills simulation and structured resident eval¬ 
uation have been implemented by residency programs [49, 50, 51, 52], A study 
by Egol et al. used a written examination known as the objective structured as¬ 
sessment of technical skill (OSATS) to assess the improvement of residents after 
a hands-on skills course and found that it was a valuable training source [51]. 
This assessment is characterized by having an experienced surgeon evaluate a 
task and rate the participant based on whether the task was: A. done correctly, 
B. partially done, or C. not done/done incorrectly A recent study by Hohn 
et al. qualitatively and quantitatively evaluated outcomes of a surgical drilling 
course through measures evaluating the ability to hit an intended target[52], A 
more sophisticated assessment of surgical skill has been reported by Karam et 
al. [49] This group implemented a one month surgical skills training experience 
for year-1 residents in their orthopaedic program. Assessment was derived by a 
combined qualitative evaluation through OSATS and quantitative evaluation us¬ 
ing a motion capture system. Using high resolution cameras and retroreflective 
markers placed on the hands of the surgeon, the researchers captured data on 
the number of discrete hand motions used to perform the desired task and the 
distance that each hand traveled to perform said task. 

In the context of orthopaedic drilling, successful task performance relates 
to creating a hole at the correct location with the correct orientation without the 
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use of excessive force, heating, over-penetration, toggle, or skiving with the drill 
to maximize efficiency and prevent harm. [40, 53, 54] The interdependency of 
task-related parameters is significant. For instance, the relationship between drill 
kinematics to force and heating are governed by drill bit type, drill RPM, bone 
parameters, feed rate, etc[55,56,57], Building off work done by Karam et al., the 
primary goal of the work presented in this chapter is to quantitatively assess or¬ 
thopaedic drilling performance. This assessment was conducted at an American 
Board of Orthopaedic Surgery-approved surgical skills training course. 


2.2 Methods 

2.2.1 Human Subjects 

This study was approved by the institutional review boards at both NMT 
(New Mexico Institute of Mining and Technology) (IRB #2014-10-004) and UNM 
(University of New Mexico) (HRRC #15-087) universities . Subjects were re¬ 
cruited during the 2015 Surgical Skills Training Course held by the Southwest 
Orthopaedic Trauma Association (SWOTA) for Year 1 orthopaedic residents in 
Phoenix, AZ. This is an ABOS residency surgical skills training course with par¬ 
ticipation from orthopaedic surgery residents from New Mexico, Arizona, Col¬ 
orado, Texas, and Utah. In all, 22 residents and 7 attending surgeons were in¬ 
cluded (9 female, 20 male, with ages ranging from 36 to 62 years, height range 
from 5 / l ,/ to 6'2" feet). 


2.2.2 Task design 

Prior to the start of the course, we asked each resident to perform a sin¬ 
gle task: drill holes through both cortices of a mechanically equivalent synthetic 
bone following a pre-defined trajectory. Following the fourth day of the training 
period, we asked each participant to repeat the task. In addition, surveys were 
administered before and after the surgical skills course. In order to character¬ 
ize performance parameters and their relation to clinic ally-relevant outcomes, a 
wide variety of parameters (as defined in Table 2.1 and Figure 2.1) were collected. 
The values obtained are listed in Tables 2.2 and 2.3. 


8 



Table 2.1: Summary of quantified measured parameters 



Parameter 

Notes 

1 

Hole Loc 

Which of four locations drilled on bone 

2 

Lj 


3 

u 


4 

l 2 


5 

Subj Num 


6 

Years 

Year in residency, where applicable 

7 

Sex 


8 

Height 


9 

Holes 

Estimated total number of holes drilled in bone 

10 

Days Since 

Days since most last bone drilling 

11 

Program Q1 

Does residency program have cadeveric training 

12 

Program Q2 

Does residency program have dedicated skills training 
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PVC 

Have they practiced with PVC bone analogs 
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OR Exp 

Have they extensive OR experience 
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Ave RPM 
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Over-penetration distance 
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25 

Ave Yaw - 0y 
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Table 2.2: Mean and standard deviation of values measured pre-test. 
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Figure 2.1: Definitions of the measured parameters in the post-test. 


2.2.3 Experimental Setup 

The parameters described in tables 2.2 and 2.3 were measured using a 
combination of motion-tracking, force, and acceleration measurements, and mea¬ 
surements taken subsequently from the bones used in the study Figure 2.2. 
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Figure 2.2: Shown here is the experimental setup used in the assessment trials. 
A surgical drill is instrumented with optical markers, accelerometer, and RPM 
sensor. The synthetic bone is attached to a fixture which is also used to measure 
applied force and to which an additional accelerometer is attached. Additional 
motion capture markers are applied to the participant to track arm and, indirectly, 
torso movement. Finally, the conventions for drill orientation (roll, pitch and 
yaw) are defined with respect to the drill body. 


A Stryker System 5 rotary handpiece was adapted with sensors to mea¬ 
sure acceleration (PCB-352C03) and rotational speed (reed switch) and motion¬ 
tracking markers for tool position and orientation. The vibration of the drill was 
logged via a datalogger, (Figures 2.3 and 2.4). 



(a) (b) (c) (d) 

Figure 2.3: (a) Stryker system 5 drill, (b) PCB-352C03 accelerometer, (c) Reed 
switch, (d) Motion-tracking system's markers. 
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A fixture was designed to hold the bone specimen (Figure 2.2. It contained 
another accelerometer of the same type and a force sensor to measure load gen¬ 
erated during drilling (Futex LSB302). Data from these sensors was measured 
using a modular DAQ system (NI cDAQ) at a sampling frequency of 10 kFIz. In 
addition to this, the arms of subjects were fit with motion-tracking markers to 
monitor arm kinematics. In particular, marker sets were oriented on the upper 
and lower arm using the greater tuberosity, lateral epicondyle, and ulnar styloid 
process as anatomical landmarks. The locations of all markers were measured 
using two Optitrack motion capture systems (one setup used eight Flex 13 cam¬ 
eras, the other used six Prime 13 cameras; two stations) with a capture frequency 
of 20Hz. 



(a) (b) (c) (d) 


Figure 2.4: (a) Futek load sensor, (b) NIcDAQ data logger, (c) Flex 13 Optitrack 
camera, (d) high current usb hub. 


Pre- and post-surveys were also administered to allow correlation between 
personal characteristics and training to measured drilling performance. Each 
subject was asked to drill five holes on the marked locations on a single bone 
mounted to a custom fixture as shown in Figure 2.2 (right). The marked locations 
produced drill trajectories of approximately 0 (normal to bone surface), +/- 25, 
and +/- 45. 


2.2.4 Data post processing 

The force, acceleration, and reed switch data was post-processed with low- 
pass filtering, as necessary, before mean, max, and standard deviation values 
were determined. Motive Body software provided with the Optitrack system 
recorded the location and orientation of the drill (6p, dp, 6y ) and the rigid-body 
kinematics of the arms. The arm rigid bodies were used to calculate elbow and 
shoulder angles. The elbow center was determined by making a best fit of the 
upper- and lower-arm vectors obtained from the respective rigid bodies. The 
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shoulder centers were then determined by moving in the known direction a dis¬ 
tance determined by Dempsters anthropomorphic table. Subsequently, trunk ro¬ 
tation relative to the orientation of the bone axis was determined. Finally, the 
authors developed forward-kinematic equations to relate the drill-body location 
and orientation to drill-tip location. The force and drill trajectories were time cor¬ 
related to the motion-tracking data and manually clipped, using force data as the 
key indicator, Figures 2.5 and 2.6. 



(a) (b) 


Figure 2.5: To find the upper-body position first, (a) motive recorded the location 
and orientation of the upper and lower arm. (b) In the post-processing proce¬ 
dure, with least square method, the location of the elbows and shoulders were 
estimated. 



Figure 2.6: The drill-tip's path calculated from the location and orientation of 
the drill with forward kinematic (blue path). The low-pass-filtered path of the 
drill-bit's location(red path) 
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2.2.5 Post hoc analysis of bone specimens 

The integrity and accuracy of drilled holes was measured directly in the 
laboratory after the conclusion of SWOTA. Digital calipers were used to mea¬ 
sure the accuracy (L^ in Tables 2.2 and 2.3) of the holes, as well as the major 
(Li) and minor (L 2 ) axis lengths. Figure 2.1. The pull-out strength was measured 
by installing 3.5 mm cancellous bone screws into the mechanically equivalent 
synthetic bones (Sawbones, Pacific Research Laboratories, Vashon Island, WA; 
Figure 2.7). Specimens were mounted to a modified angle vise that allowed po¬ 
sitioning the screws in line with the test actuator. Figure 2.8. The bones were 
cyclically loaded using an MTS 858 Mini Bionix II loading system using a proto¬ 
col described in Ezechieli et al. [58] Briefly, the protocol is as follows: 

1. Screw pre-tensioned with a constant load of 60 N for 60 s 

2. Adding 30 s of dwell 

3. Followed by 500 cycles between 60 N and 250 N at 0.75 FIz 

4. Force was ramped down to zero 

5. Followed by a 30 s break 

6. A pull-out test was subsequently carried out at 1 mm/s to failure 



Figure 2.7: The pullout test setup to measure the maximum strength of the stan¬ 
dard screws inside each hole. 
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Figure 2.8: Various types of failure of the bone in the pullout test 


Figure 2.9 shows the general flow of the data in one chart. The data from 
each part of the setup, marked with a time stamp from the CPU, and then the 
time stamp used to synchronize and clip all of the data. 
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Figure 2.9: The general flow chart of data from various channels and devices, 
stamped with the CPU time for each device. 
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2.3 Results 


The assessment data were obtained through a combination of post-processing 
data obtained by the motion-tracking system, force sensor, accelerometers, reed 
switch, and subsequent measurements made on the bones used in the study, as 
summarized in Tables 2.2 and 2.3. To assess the interdependency of the mea¬ 
sured parameters, a covariance matrix for all variables was generated (Figure 
2.10). Non-trivial (e.g. sex with height) correlations whose magnitude exceeded 
0.5 were (correlation values listed taken from post data set): 

• Maximum penetration distance with force (SD, 0.72) 

• Maximum penetration distance with force (max, 0.73) 

• Maximum penetration distance with force (mean, 0.91) 

• Drill toggle with drill roll angle (mean, 0.70) 

• Drill RPM with force (mean, 0.53) 




Figure 2.10: A visualization of how various parameters were found to co-vary 
is given by color-mapping correction coefficients. Correlations on color map 
toword blue include negative correlations. Axes are defined by variable in ta¬ 
ble 3.1. Plot (a) is the result from the test before the SOWTA training and plot (b) 
is the result from the test after the SOWTA training. 


To have a better understanding of the data, the result of each parameter 
plotted separately respect to the user id. The best performance for the majority 
of parameters belonged to experts. Figure 2.11. 
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(a) (b) (c) (d) 

Figure 2.11: Plot of the ranked data for average (a) force, (b) mean roll of the drill, 
(c) maximum penetration of the drill and, (d) final minor axis of the drilled hole 
vs the hole number. For each plot we checked the id number of the winner and 
see if it was an expert or a novice. 


2.3.1 Pre and Post Survey 

Full summery results are included in the Appendix A. Key parameters are 
included in Tables 2.2 and 2.3. The participant's feedbacks will be consider in the 
future studies. 


2.4 Conclusions 

This study used various sensors to measure the parameters addressed by 
the trainers as the most critical factors in bone drilling. To find the effects that 
these parameters have on each other, the measured parameters were expressed 
as scalars and their covariances were calculated. Some immediately evident re¬ 
sults include that the covariance matrices are not drastically different before and 
after the training, nor are the parameters themselves as can be seen by compar¬ 
ing Tables 2.2 and 2.3. This may not be a surprising result given that the training 
which occurred over the intervening time covered many skills, with no particular 
emphasis on drilling. The data suggests that when taken as a group, the experi¬ 
enced surgeons tested here did not, on average, outperform the residents collec¬ 
tively. Flowever, when looking at clinically relevant parameters independently, 
an interesting pattern does seem to provide some insight into the effects of addi¬ 
tional experience. The very best performers out of the total 29 participants, when 
looking at one performance aspect at a time, came from the experienced surgeon 
subset with respect to: maintaining steady force (highest minimum before exit¬ 
ing the distal cortex), minimizing over-penetration, minimizing both the major 
and minor axis diameters, minimizing toggle and drill vibration. Additionally, 
an experience surgeon was nearly (second) best in terms of accurately targeting 
the mark on the distal cortex and (third) best in terms of the force before failure 
during subsequent screw pull-out testing. 
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No parameters related to arm or torso kinematics were found to vary in a 
patterned manner with respect to other parameters. This may be due to the vary¬ 
ing height and dominate handedness of participants as well as the observation 
that surgeons were quite unique in their chosen stance. 

Limitations of this study include the fact that three residents had to leave 
the event early and could not be retested, that the pool of experienced surgeons 
was small and perhaps not representative (some had moved into specializations 
that did not involve significant amounts of bone drilling), and that a significant 
fraction of the motion-tracking data had to be excluded due to occlusions which 
resulted from the manner in which surgeons chose to position themselves dur¬ 
ing the task. Additionally, many surgeons commented in the post-survey that 
either the bones felt too hard or that the drill bit felt dull. This is perhaps not 
surprising given that the synthetic bones used in this trial were composite bone 
models, rather than the compressed foam-type that is more commonly used for 
training. Planned future work includes major revisions to the instrumentation 
to reduce post-processing time and expanding the study to include a larger par¬ 
ticipant pool. We expect that this will further elucidate the differences between 
novice and experienced orthopaedic surgeons. Ultimately, we hope to translate 
the outcomes into improved surgical training methodologies. 
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CHAPTER 3 


DESIGN OF AN AUGMENTED HAPTICS DEVICE FOR 
ORTHOPAEDIC SURGICAL SKILLS TRAINING 


3.1 Introduction 

Orthopaedic surgery requires a wide range of technical skills that are mainly 
acquired through an apprenticeship model. Recently the American Board of Or¬ 
thopaedic Surgery agreed to increase skills training requirements for orthopaedic 
surgical residents [59]. Specific skills can be practiced on cadavers, synthetic ma¬ 
terials, and animals. However, these materials can be prohibitively expensive 
[60]. Furthermore, none of these training methods constitutes a perfect analog 
to live surgery on humans. Even the simple task of drilling a hole in the bone 
prior to fixation requires that the surgeon simultaneously: (1) avoid skivvying 
(slipping) off the bone when beginning a hole, especially when oblique, (2) main¬ 
tain a desired orientation vs. toggling, (3) drill rapidly to avoid overheating the 
bone [49], (4) anticipate the changes in the mechanical properties of the bone re¬ 
lated to its non-uniform profile [50], (5) drill through the entire width of the bone 
for a bicortical screw placement, and yet (6) avoid over-penetration. These as¬ 
pects are greatly complicated by the wide range of bone geometry and material 
properties that relate to age, sex, and disease state.The focus of this work is train¬ 
ing orthopaedic residents in the placement of bicortical screws. These screws are 
meant to span the entire width of the bone, engaging on both cortices. 

Significant recent work has explored the use of virtual environments (VEs) 
and robotic tools to simulate various surgical tasks, including those relevant to 
orthopaedic surgery [51, 40, 52, 61, 54, 23, 62, 47], For those VEs where a manip¬ 
ulation task is involved, haptic elements are often incorporated. One of the first 
VEs for orthopedic surgery was developed by Dar Tsai et al. [55]. These simu¬ 
lations have been generally received with cautious optimism, despite significant 
limitations that often exist in fidelity [52], For instance Fang et al. received feed¬ 
back suggesting that their VE did not allow surgeons to know the strength of the 
drilling bone or the distance between the drill-tip and the target location on the 
temporal bone [ 57 ]. In fact, studies have shown that axial forces range widely 
from task to task, but can easily exceed 200 N [63]. Moreover, once the drill has 
penetrated some distance into the bone, torque feedback is evident and critical 
to maintaining a desired trajectory, estimating hole depth, and monitoring bone 
properties. As such, those VEs that have previously relied on haptic feedback are 
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intentionally limited to procedures that required low force and/or no torque [52], 

Augmented Environments (AEs) have been used for a wide range of pur¬ 
poses [58, 49, 58], including related tasks such as the feeling of skin when insert¬ 
ing a needle into the body [51]. Likewise for bone drilling, the combined use 
of a physical synthetic bone for realistic torque and axial force magnitudes with 
additional force from a custom haptic interface seems promising. Critically, the 
ability to arbitrarily affect net force enables the simulation of anatomic variability 
of bone. Moreover, if this method is shown to perform well, it might be extended 
to allow the use of low-cost (e.g. PVC pipe) materials, which are currently used 
in locations where funding is particularly limited. 

A large fraction of the previous research in augmented haptics utilizes 
real and virtual elements dealt with separately. For instance, forces on one hand 
would come from interaction with a real object while forces on other hand would 
come from interaction with a virtual environment by means of haptic technology. 
Here we focus on applications where a the task is best suited to physical coinci¬ 
dence of virtual and real aspects - both in terms of contact/force and visual ele¬ 
ments. As one example, Kurita et al. combined an elastic gel with force feedback 
to render visco-elastic behavior [64], Similarly, Kuchenbecker et al. superposed 
high-frequency force transients during impact with physical objects [65]. Jeon 
et al. made an augmented-haptic model for breast cancer palpation, combining 
a silicone breast model with a virtual tumor [66]. In another application, a vir¬ 
tual sphere was rendered coincident to a rubber ball in order to compensate for 
actuator saturation during deep penetration [32]. 

The specific task of modeling cutting/drilling bone has been pursued pre¬ 
viously, including first-principles approaches [52, 63]. For a given drill bit, and 
where rotational speed material properties are constant, the model in [52] pre¬ 
dicts an axial force that scales linearly with feedrate. Yanping et al. used the 
damping model for simulating bone sawing [40]. Tsai et al. used an exclusive 
model for the haptic interaction orthopedic simulator built upon the feed rate 
and torque of the drill [66]. Lee et al. modeled the drilling of the bone with re¬ 
spect to the rpm and feed rate [50]. In some studies force is modeled linearly with 
respect to feed rate [67], 


3.2 Methods 

This study attempts to combine positive aspects of synthetic bones with 
virtual reality as an augmented reality environment. Residents will have different 
perception from the applied force feedback from the haptic environment while 
drilling in the same bone. This difference in feeling trains residents to drill in 
various bones with different densities. Therefore the rendered force is based on 
the difference of the force in the various cross-sections of synthetic bones. 
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3.2.1 Hardware Design and Calibration 

To characterize the mechanics of drilling into long bone, a system of simul¬ 
taneously recording kinematics and force was required. A Phantom Premium 
(high force model. Geomagic, Rock Hill, SC) already capable of measuring gim- 
bal location, was modified to connect to a surgical drill while capturing pitch and 
yaw angles (Figure 3.1). 



Figure 3.1: (a) Phantom Premium 1.5HF mounted to the Stryker System 5 drill 
via a customized gimbal. Yaw and pitch angle are captured using an Arduino 
Mega and two Avago HEDM-5500 encoders, (b) The forward-kinematics for the 
modified device use the length and angle definitions shown. The joint angles 
{q\,q 2 , • • •, ^ 5 } of device are defined with the device in the home configuration. 


In order to find the absolute location and orientation of the drill, the forward- 
kinematics needed to be recalculated taking into account changes to the structure 
of the haptic device. We define relative joint angles {q 1 , qi,... ,q$) as shown in 
Figure 3.1(b) in order to capture the kinematics associated with the standard prox¬ 
imal structure as well as the distal structure including the custom gimbal and 
drill geometry (Figure3.2). Accordingly, the drill-tip position and orientation can 
be determined using the homogeneous transformation matrix (using notation by 
Spong [ 68 ]): 


H = Rof Ml Rot zm Trans X/Ll Rot ZrCj3 ... 

■ ■ ■ Trans XiL2 Rot X/Cj4 Rot ZrCj5 Trans X/L3 


(3.1) 


Because of the relatively stiff surface of the bone and the desire to apply no 
force before contact, forward-kinematics must be calculated with accuracy and 
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precision. An external metrology system (an 8-camera Optitrak Flexl3 system) 
was used to verify drill-tip location in various configurations. No errors exceed¬ 
ing 1 mm were observed in any configuration. Note that this is not much greater 
than the uncertainty introduced by the metrology system itself (0.6 mm in plane, 
0.7 mm out of plane). 



Figure 3.2: (a) 3D model of the gimbal. (b) The schematic of the drill 


The forward kinematic model was then used to derive the manipulator 
Jacobian. Its terms (partial derivatives) were obtained using Maple and are in¬ 
cluded in Appendix C.l. To verify the accuracy of the result, the actual versus 
commanded forces were compared at various configurations and found to agree 
well (Fig. 3.3). 


24 










Figure 3.3: To validate the Jacobian matrix, the output force of the device was 
measured by a Shimpo force gage, relative to the commanded value at various 
configurations. 


3.2.2 Bone Drilling Mechanical Model 



(a) (b) (c) (d) 


Figure 3.4: Cross-sections of two types of synthetic bone were used for the pur¬ 
pose of modeling each with differing thickness of cortical and cancellous (trabec¬ 
ular) bone, (a) Hard bone with cancellous , (b) hard bone without cancellous , (c) 
softbone with cancellous and, (d) softbone without cancellous 


Two major parameters related to the mechanics of drilling are thickness 
and hardness of the cortical layer. Synthetic bone specimens with these variations 
were prepared (Fig. 3.4). Multiple holes were drilled in each specimen using the 
apparatus shown in Figure 3.5. The force during drilling was measured with 
a Futek force sensor (model LLB400) and a NI data logger (cDAQ-9184) with a 
sampling frequency of 100 Hz. The Phantom Premium device (for {cj i ... cjj }) 
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in conjunction with an Arduino Mega (for ^4 and ^ 5 ) were used to measure the 
location of the drill bit with a frequency of 50 Hz. The kinematic data was post- 
processed using a 6 th order Butterworth filter with 0.1 Hz cut-off frequency. 



Figure 3.5: Test setup for measuring the force during the drilling. 


Initially, a linear, velocity-dependent model was used to fit the data, as per 
the recommendations in the literature mentioned previously. This effort resulted 
in a poor fit to the data. Instead a phenomenological approach was taken, using 
the non-linear objective function: 

A O /I O 

f = ( 14 .x -\- CI 3 X T-... T" ciq T" fqx -\- b^x -F ... T~ (3-2) 


with the corresponding multiple regression model: 
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(3.3) 


where n is the number of data points gathered in all trials for a particular bone 
specimen and x is defined normal to the surface of the bone with the origin at 
the surface and with up (away from bone) being positive. In order to create sim¬ 
ulate the mechanics of drilling into bone with sufficient accuracy for the sake 
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of training, we sought to create a phenomenlogical model with optimized accu¬ 
racy, with focus on position- and velocity-dependent parameters. Position and 
velocity have been shown previously to critically underlie drilling mechanics 
[11]. Favoring simplicity, in light of the need for real-time performance, and it¬ 
erative algorithm started with a first order polynomial model in x and x, (F = 
a\x + flo + b\X + bo). Higher order terms of both x and x were added and a re¬ 
gression was run to calculate the root-mean-square error (RMSE). This RMSE was 
compared to that from the previous step (simpler model) until the improvement 
in error dropped to below 10%." The final models produced were found describe 
the respective bone types with the RMSE errors (given in Figures 3.6 - 3.9). 


3.3 Results 
3.3.1 Force Modeling 

Using equation(3.3), model parameters were obtained for each bone sam¬ 
ple (Table 3.1). As expected, the various trials for the same sample showed 
some variability, but were well described by the respective models nonetheless 
(Figs. 3.6, 3.7, 3.8, 3.9). By contrast, the models (surfaces) found for the various 
samples were quite distinct from each other (Fig. 3.10), suggesting that the sam¬ 
ples tested are indeed mechanically distinct. It can be further conjectured that the 
role of the thickness of the cortical and cancellous regions play a significant role. 

It should be noted that the domain used in fitting the surface is not thor¬ 
oughly explored by the data obtained. The extra domain is somewhat unavoid¬ 
able as certain combinations are difficult to attain, such as high velocity at the 
initial position. This is not a significant limitation, as such a combination is not 
relevant to this application. 
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parameter 

Sample A 

Sample B 

Sample C 

Sample D 


-9010 

-5220 

-5080 

-19060 

a 3 

-11900 

-6280 

-7030 

-17100 

Cl 2 

-5440 

-2730 

-3080 

-5200 

U\ 

-939 

-516 

-454 

-601 

Uq 

0.351 

10.9 

2.10 

3.55 

h 

792 

-14800 

-1880 

-48900 

h 

393 

-19100 

-1930 

-24700 

b 2 

-42.8 

4330 

-583 

-3460 

h 

-28.8 

-95.4 

-55.9 

-130 

ao 

0.351 

10.9 

2.10 

3.55 


Table 3.1: Model (surface) parameters, where the units of each term are such 
that when multiplied by the corresponding position-, velocity-, or offset-based 
parameter (e.g. x) the result is in Newtons. Note that the kinematic equations 
were derived using inches as the length unit to match design and manufacturing 
around imperial units. 
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Figure 3.6: Plot of the force model with respect to the measured data sets (mul¬ 
tiple trials) are shown for the sample comprised of hard bone with a cancellous 
region (Sample A), RMSE 8%. The model is shown to fit the data well; how¬ 
ever, it is critical to note that the model extends well beyond the domain actually 
sampled. This is also true for the models shown in the immediately following 
figures. 
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Figure 3.7: Plot of the force model with respect to the measured data sets (multi¬ 
ple trials) are shown for the sample comprised of hard bone without a cancellous 
region (Sample B), RMSE 9%. 



Figure 3.8: Plot of the force model with respect to the measured data sets (mul¬ 
tiple trials) are shown for the sample comprised of soft bone with a cancellous 
region(Sample C), RMSE 11%. 
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Figure 3.9: Plot of the force model with respect to the measured data sets (multi¬ 
ple trials) are shown for the sample comprised of soft bone without a cancellous 
region (Sample D), RMSE 13%. 
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Figure 3.10: Comparison of the four force models, each resulting from one of the 
bone samples tested, also shown in figures 3.6 - 3.9. 


Figures. 3.11, 3.12, 3.13 and 3.14 compare multiple trials for each material. 
Although the error in some areas is considerable, most of these areas are not 
reachable for location and velocity. 
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Figure 3.11: Comparing the four force models with respect to the measured data 
sets are shown for the sample comprised of hard bone with a cancellous region 
(Sample A) 



Figure 3.12: Comparing the four force models with respect to the measured data 
sets are shown for the sample comprised of hard bone without a cancellous region 
(Sample B) 
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Figure 3.13: Comparing the four force models with respect to the measured data 
sets are shown for the sample comprised of soft bone with a cancellous region 
(Sample C) 



o 


Figure 3.14: Comparing the four force models with respect to the measured data 
sets are shown for the sample comprised of soft bone without a cancellous region 
(Sample D) 
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Figure 3.15: Plot of the average force model with respect to the measured data 
sets (all trials) are shown for the sample comprised of hard bone with a cancellous 
region (Sample A) 



Figure 3.16: Plot of the average force model with respect to the measured data sets 
(all trials) are shown for the sample comprised of hard bone without a cancellous 
region (Sample B) 
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Figure 3.17: Plot of the average force model with respect to the measured data 
sets (all trials) are shown for the sample comprised of soft bone with a cancellous 
region (Sample C) 



Figure 3.18: Plot of the average force model with respect to the measured data sets 
(all trials) are shown for the sample comprised of soft bone without a cancellous 
region (Sample D) 


35 





Figure 3.19: The modified haptic device (right) is shown with the force-gage- 
embedded bone fixture (bottom) and data acquisition system (right). In order to 
obtain the kinematic parameters ( x, x), the raw enocder values were read using a 
custom C++ script and then fed into the inverse kinematic equations we derived. 
The force sensor was driven by and read using an National Instruments data 
acquisition system. 
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3.3.2 Test Evaluation 


To check the performance of the method, a custom chassis was constructed. 
As a first step to check the modified Jacobian, a vertically sliding table was de¬ 
signed and the robot set on top of it. The force sensor under the sliding table, 
measured the vertical load that was applied from the robot during the drilling 
(Figs. 3.20). 



Figure 3.20: The test setup to check the modified Jacobian. By putting the robot 
on a force sensor in a vertically slider table, during the evaluation test, the ap¬ 
plied force from the robot was continuously measured and compared with the 
calculated force. 


As an experiment, in the case where bone B is meant to feel like A, the 
algorithm would: 

• Measure the location x, and the speed x. 

• Estimate the current force associated with drilling into B, given the current 
kinematics. 

• Calculate the necessary difference in force to make B feel like A. 

• Apply this force to the user. 
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Figure 3.21 shows the result of the validation proses. 



Figure 3.21: The applied force in validation study. 


To make the system stable, the reference force trajectory was overlaid on a 
signal ramp function. 


3.4 Conclusions 

The force model which was used in this study is a function of drill bit's 
position and velocity relative to the bone surface. This model is more complex 
than those presented by others in the prior work described and better fit the data 
we obtained, by a significant margin, than these simpler models. It may be that 
the position dependence relates to frictional forces, which increase with depth. In 
particular, these simpler models did not contain the non-linear terms present in 
our model. This suggests that better fit provided by our model points to drilling 
mechanics being more complicated than previously thought. It should again be 
mentioned that the cross-sections are non-uniform. 

These results suggest that accurate predictive models of drilling can be 
obtained. These models may be used to augmented one bone sample, to "feel" 
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like some target bone type. To effectively harden the cortical or cancellous re¬ 
gion, this model can predict the necessary increase in force against the drilling 
direction, and visa versa. 

Future work will extend the types of bone tested to include flat, irregu¬ 
lar, and short bones as well as cadaveric samples. Future work will investigate 
the quantitative and perceived realism of augmented bone drilling through user 
studies with surgical residents and experts. This preliminary success does, in 
fact, suggest that low-cost (e.g. PVC pipe) materials might be used to reduce cost 
without sacrificing fidelity. 

Future work might also explore the use of augmenting torque feedback 
during drilling. In addition, it may prove to be useful to intentionally distort 
what is felt to exaggerate certain effects for the purpose of training (e.g. making 
the bone surface especially slippery). Force perturbation might elicit reductions 
in toggling, decreased surface friction may improve tactics related to avoiding 
skivvying, and so forth. We will also explore rendering soft tissue that would be 
adjacent to the bones, further increasing the realism of the simulation. 
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CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 


4.1 Summary 

This thesis explored the potential benefits of using robotics tools to im¬ 
prove training methods in surgical drilling. Modern tools such as motion cap¬ 
ture systems and haptic robots may enhance the precision and fidelity of training 
tools. Even subtle improvements, might play a significant role in the utility of 
bone drilling training. Moreover, improvement to training may significantly re¬ 
duce risk during surgical procedures. Here, through quantification of mechanical 
and "human moto" parameters we present an improved understanding of how 
a surgeons control of a drill affects relevant outcomes. 


4.1.1 Skill Assessment 

Clearly a surgeons ability to correctly execute a surgical path relies on 
some combination of dexterity, clinical knowledge, and prior experience. It is 
worth considering to what extent skill assessments such as those presented here 
specifically address each of these areas. As there was no significant improve¬ 
ment in recorded parameters during the traditional training, this study suggests 
adding more effective and modern methods. 

The fact that the performance parameters measured were not significantly 
different over the course of the focused training course like SWOTA, suggests 
that either substantial improvements in the training method is required or the 
assessment tool must be revised to better capture the differences between and 
novice and expert performance. 


4.1.2 Augmented Haptics 

The use of virtual reality for surgical training could yield many exciting 
benefits. However, limitations exist in the accuracy of the mechanical models as 
well as the capability of existing commercial haptic devices for surgical training. 
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To improve realism, VR systems for bone drilling need a more accurate mechan¬ 
ical model. Here we were able to demonstrate an accurate mechanical model. 
The use of augmented haptics was shown to critically extend the capability of 
standard, impedance-type haptic devices for tasks like training bone drilling. 


4.2 Future Work 
4.2.1 Skill Assessment 

In this thesis a case for additional studies on the skill assessment of bone 
drilling was presented. This study demonstrated that there are limitations in the 
sensitivity of the measurement tools used here relative to the surgical skill eval¬ 
uation and/or specificity of training events like SWOTA. It stands to reason that 
focused training on a particular skill, such as avoiding over-penetration, would 
result in more significant/rapid improvements in that particular aspect. Clearly 
this is difficult to accomplish at these events given the wide range of skills that 
need to be acquired in the relatively short time. Some of the assumptions, such as 
the effect of vibration and the angle of penetration, did not show a strong correla¬ 
tion with over-penetration or other clinically important parameters. The current 
trial design is too bulky, but before being able to make the design simpler and 
more compact, we must understand which parameters may be neglected. It is 
necessary to conduct more tests and prioritize the measured parameters. More 
than simplification, having more data and participants can significantly improve 
the quality and accuracy of the data. 

Having clear definitions for parameters can help to generate a more accu¬ 
rate index. This study used traditional definitions for the measured parameters; 
as an example, we assigned the "expert surgeon" class to a practitioner past their 
residency. Titerature showed that spending more time in the OR does not neces¬ 
sarily make an expert surgeon, therefore making better definitions of parameters 
can improve the quality of the training significantly. Experts participated in the 
studies during different conferences and events, and occasionally the team had 
difficulties explaining the task clearly due to the lack of time subject focus. Not 
being consistent with the explanation of the task, as well as letting participants 
communicate during the trial may have added error. Defining a more routine 
task may magnify the difference between experts and novices, and experiment¬ 
ing with bones of various mechanical properties is also recommended. 

Future work might include the use of complex Neural Network to further 
elucidate the connection between clinically relevant outcome of performance pa¬ 
rameters. So far only a simple, single hidden layer Neural Network was used but 
this study did not produce any significant results suggesting that a more compli¬ 
cated Neural Network maybe required. Data-mining and big-data methods can 
also be beneficial for future work[69]. 
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4.2.2 Augmented Haptics 

The physical model can easily move compared to the registered location 
in the VR coordinate system. Therefore, adding a physical collision detection 
system can help with the performance of the device, providing superior accuracy 
over a location estimate. As another approach, adding an image processing base 
setup can continuously register the location of the model via virtual coordinates. 

The force model presented here for drilling has only modest error. Addi¬ 
tional work related to fine tunning the model can result significant improvement. 
Additionally conducting trials on wider range of materials and drill bit could 
expand the utility of this model. 

To generate a more realistic virtual reality it is necessary to add a second 
and third dimension to the model. Furthermore, adding torque feedback maybe 
critical to realism, although it is worth noting that most commercial haptic de¬ 
vices are incapable of providing the level of torque required, if any at all. The 
available haptic devices capable of torque feedback tend to have much smaller 
work space than their equivalent without torque feedback. 
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APPENDIX A 


SURVEYING RESULT FROM NOVICE PARTICIPANTS 
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Id Record Name Date Created Date Updated Record Stat 1. Years in 12. Sex 3. Handedn 


1 Record_l 

5/14/15 9:06 

5/14/15 13:55 Saved 

5 

1 

1 

2 Record_2 

5/14/15 10:57 

5/14/15 13:55 Saved 

1 

1 

1 

3 Record_3 

5/14/15 11:15 

5/14/15 11:19 Saved 

1 

1 

1 

4 Record_4 

5/14/15 11:33 

5/14/15 11:54 Saved 

1 

1 

1 

5 Record_5 

5/14/15 11:55 

5/14/15 11:58 Saved 

1 

0 

1 

6 Record_6 

5/14/15 12:09 

5/14/15 12:13 Saved 

1 

1 

1 

7 Record_7 

5/14/15 12:13 

5/14/15 12:36 Saved 

1 

1 

1 

8 Record_8 

5/14/15 12:42 

5/14/15 12:45 Saved 

1 

1 

1 

9 Record_9 

5/14/15 12:45 

5/14/15 12:48 Saved 

1 

1 

1 

10 Record_10 

5/14/15 13:49 

5/14/15 14:26 Saved 

1 

1 

1 

11 Record_ll 

5/14/15 13:58 

5/14/15 14:27 Saved 

1 

0 

1 

12 Record_12 

5/14/15 14:43 

5/14/15 14:48 Saved 

1 

0 

1 

13 Record_13 

5/14/15 14:49 

5/14/15 14:53 Saved 

1 

1 

1 

14 Record_14 

5/14/15 15:04 

5/14/15 15:07 Saved 

1 

0 

1 

15 Record_15 

5/14/15 15:09 

5/14/15 15:12 Saved 

1 

1 

1 

16 Record_16 

5/14/15 15:22 

5/14/15 15:25 Saved 

1 

1 

1 

17 Record_17 

5/14/15 15:25 

5/14/15 15:29 Saved 

1 

1 

0 

18 Record_18 

5/14/15 15:37 

5/14/15 17:01 Saved 

1 

1 

1 

19 Record_19 

5/14/15 15:58 

5/14/15 17:01 Saved 

1 

0 

1 

20 Record_20 

5/14/15 16:02 

5/14/15 17:18 Saved 

1 

1 

1 

21 Record_21 

5/14/15 16:20 

5/14/15 17:17 Saved 

1 

0 

1 

22 Record_22 

5/14/15 16:53 

5/14/15 17:17 Saved 

1 

1 

1 

23 Record_23 

5/14/15 17:14 

5/14/15 17:17 Saved 

1 

1 

1 
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4. Height (in) 


5. How many holt 6. How mar 7. How mar 8. Does you 9. Do you b If 1, in whal 10. Does yo 


72 

70 

100 

1 

10 

1 

1 

90 

0 

1 

1 

1 nan 

1 

0 

73 

5 

0 

2 

1 

1 1,2 


1 

71 

3 

1 

28 

0 

1 

1 

0 

68 

10 

10 

3 

1 

1 

1 

0 

70 

230 

1 

30 

0 

1 1,2 


0 

75 

10 

1 

4 

1 

1 1,2 


1 

70 

2 

2 

5 

0 

1 1,2 


0 

75 

10 

100 

2 

1 

0 nan 


1 

66 

2 

0 

5 

1 

1 1,2 


1 

67 

1 

0 

2 

1 

1 

1 

1 

68 

1 

1 

4 

0 

1 

1 

0 

69 

3 

0 

1 

1 

1 

1 

1 

64 

10 

0 

7 

1 

1 

1 

1 

69 

4 

0 

14 

0 

0 nan 


1 

68 

20 

5 

1 

1 

1 

1 

1 

73 

10 

0 

1 

1 

0 nan 


1 

70 

5 

0 

10 

1 

1 

1 

1 

69 

10 

10 

2 

0 

1 

1 

1 

69 

20 

0 

2 

1 

1 

1 

1 

64 

20 

0 

2 

1 

1 1,2 


1 

72 

1 

0 

90 

0 

1 1,2,3,4,5 


0 

72 

1 

0 

90 

0 

1 

1 

0 
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11. Does yc 12. Does yo 13. Does yc 14. How are your sui other: 
nan nan nan nan nan 

1 0 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

111 c. Objective Structured Assessment of Technical 
1 0 Id. other Ot 

1 0 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

111 d. other Month training was second, 

1 1 1 a. Faculty Evaluation 

111 d. other all three 

0 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 

1 1 Id.other a and c 

1 1 1 a. Faculty Evaluation 

1 0 1 a. Faculty Evaluation 

1 1 1 a. Faculty Evaluation 


15. Do you 16. Do you 
nan nan 


0 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

1 

0 

1 

1 

1 

1 

1 


0 

1 

0 

0 

0 

1 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

1 

1 

0 

0 
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17. Do you a. If 1, did y 18. 

Do you a. 

If 1, how 19. 

On a sc; 20. (21. Participant Number 

Addition 

1 nan 


1 nan 

8 

8 

5 

1 Ot in resi 

1 

1 

1 

400 

8 

8 

7 

2 UNM sti 

1 

1 

0 


7 

8 

7 

5 

1 nan 


1 

150 

5 

5 

3 

6 Done bo 

1 

1 

1 

20 

6 

6 

6 

7 

1 

1 

0 


6 

4 

3 

3 

1 

1 

1 

200 

7 

7 

7 

8 Multiple 

1 

1 

1 

100 

7 

7 

5 

9 

1 

1 

1 

1000 

5 

4 

2 

4 

1 

1 

1 

8 

7 

6 

8 

10 

1 

1 

0 


4 

4 

2 

11 

1 

1 

0 


5 

4 

6 

13 

1 

1 

0 


6 

5 

6 

12 

1 

1 

0 


3 

3 

4 

14 

1 

1 

1 

100 

8 

4 

4 

17 

1 

1 

0 


5 

6 

6 

15 

1 

1 

1 

1000 

5 

5 

3 

18 

1 

1 

0 


4 

6 

3 

16 

1 

1 

1 

24 

8 

8 

8 

19 

1 

1 

0 


8 

8 

6 

23 

1 

1 

0 


4 

5 

6 

20 Had resi 

1 

1 

0 


3 

3 

4 

21 

1 

1 

1 

4 

5 

6 

6 

22 
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Id 

idency 


Record Nan Date Create Date Updat Record Stat 1. Do you h a. If 1, how 2. On a seal 


8 Record_8 ######### ######### Saved 1 100 8 

4 Record_4 ######### ######### Saved 0 7 

th models First year in this residency, did one year general prior to this one E6 - 0 handed 


19 

Record_ 

.19 

######### 

######### 

Saved 

1 

40 

6 

1 

Record_ 

.1 

######### 

######### 

Saved 

0 


5 

14 

Record_ 

.14 

######### 

######### 

Saved 

1 

200 

7 

7 

Record_ 

.7 

######### 

######### 

Saved 

1 


7 

11 

Record_ 

.11 

######### 

######### 

Saved 

1 

2500 

6 

3 

Record_ 

.3 

######### 

######### 

Saved 

0 


7 

2 

Record_ 

.2 

######### 

######### 

Saved 

0 


4 

22 

Record_ 

.22 

######### 

######### 

Saved 

0 


6 

5 

Record_ 

.5 

######### 

######### 

Saved 

0 


5 

26 

Record_ 

.26 

######### 

######### 

Saved 

0 


4 

9 

Record_ 

.9 

######### 

######### 

Saved 

1 

100 

8 

16 

Record_ 

.16 

######### 

######### 

Saved 

0 


6 

27 

Record_ 

.27 

######### 

######### 

Saved 

0 


4 

13 

Record_ 

.13 

######### 

######### 

Saved 

0 


8 

10 

Record_ 

.10 

######### 

######### 

Saved 

1 

24 

8 

18 

Record_ 

18 

######### 

######### 

Saved 

0 


5 

6 

Record_ 

.6 

######### 

######### 

Saved 

0 


4 

12 

Record 

12 

######### 

######### 

Saved 

1 

4 

6 


48 



3. On a seal 4. On a seal Previous Recurrent Res Additional Ote Id 


Record Nar Date Create Date Updat 


8 

8 

2 

30 

8 Record, 

_8 

######### ######### 

7 

7 

5 

26 New rigid defii 

4 Record- 

_4 

######### ######### 

6 

6 

7 

43 

20 Record- 

-20 

######### ######### 

5 

6 

3 

24 First hole - tab 

1 Record- 

-1 

######### ######### 

7 

5 

8 

38 After c, drill bi 

15 Record- 

.15 

######### ######### 

7 

8 

9 

31 

7 Record- 

_7 

######### ######### 

7 

5 

4 

35 

11 Record- 

-11 

######### ######### 

7 

8 

10 

27 

3 Record- 

_3 

######### ######### 

4 

3 

11 

25 Problem with 

2 Record- 

_2 

######### ######### 

5 

6 

13 
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5 

4 

12 

28 

5 Record- 

_5 

######### ######### 

4 

5 

14 

50 Used step stoc 

27 Record- 

-27 

######### ######### 

7 

7 

17 

33 Right side 

9 Record- 

_9 

######### ######### 

6 

7 

15 

39 

17 Record- 

-17 

######### ######### 





14 Record- 

.14 

######### ######### 

6 

5 

16 

51 starting on B 

28 Record- 

.28 

######### ######### 

7 

8 

19 

36 Drill bit is gett 

13 Record- 

.13 

######### ######### 

8 

8 

23 

32 

10 Record- 

.10 

######### ######### 

6 

7 

20 

41 Standing on st 

19 Record- 

.19 

######### ######### 

4 

5 

21 

29 

6 Record- 

_6 

######### ######### 

6 

7 

22 

34 
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Record Stat 1. Do you fr a. If 1, whic 2. On a seal 3. On a seal 4. On a seal 5. Do you b 6. Do you f(7. Do you b 


Saved 

1 markers on 

8 

8 

8 

1 

1 

1 

Saved 

0 

8 

6 

7 

0 

1 

1 

Saved 

0 

6 

7 

6 

1 

1 

1 

Saved 

1 tracker ball 

5 

4 

6 

1 

0 

1 

Saved 

1 drill was du 

8 

9 

8 

1 

1 

1 

Saved 

0 

7 

7 

8 

1 

1 

1 

Saved 

1 location am 

5 

9 

5 

1 

1 

0 

Saved 

0 

7 

7 

8 

1 

1 

1 

Saved 

1 markers on 

4 

4 

3 

1 

1 

1 

Saved 

0 

5 

3 

5 

1 

1 

1 

Saved 

0 

4 

4 

5 

1 

1 

1 

Saved 

0 

8 

8 

7 

1 

1 

1 

Saved 

0 

6 

7 

7 

1 

1 

1 

Saved 

0 

6 

6 

4 

1 

1 

1 

Saved 

0 

4 

6 

4 

1 

1 

1 

Saved 

0 

5 

8 

8 

1 

0 

1 

Saved 

1 hand on dri 

8 

8 

8 

1 

1 

1 

Saved 

1 heigntwas 

6 

6 

7 

1 

0 

1 

Saved 

0 

5 

6 

6 

1 

1 

1 
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a. If 1, in wt 8. Do you h Sides Drillei Previous Re Current Res Additional Otes 



1 

Markers obstructing av 

2 

30 


1 

0 

5 

26 


1 

First Time 1 

7 

43 

1,2 


More clearly define obj 

3 

24 

1,2 


First Time 1 

8 

38 

1,2 



9 

31 Different side 



Wires and f First Time l 

4 

35 

1,2 



10 

27 


1 


11 

25 


1 

Getting drill more 

12 

28 

1,2,3,4,5 


First Time 1 

14 

50 

1,2,3,4,5 


First Time l 

17 

33 


1 

First Time 1 

15 

39 


1 

First Time l 

18 

37 

1,2,3,4,5 


First Time 1 

16 

51 


1 

First Time 1 

19 

36 


1 

Limited abil First Time 1 

23 

32 

1,2 


Adjustable First Time 1 

20 

41 

1,2,3,4,5 



21 

29 
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APPENDIX B 


MATLAB CODES AND FUNCTIONS 


B.l Skill assessment project 
B.1.1 Main file 1 


1 

clc 


2 

clear all 


3 

[data_lv,h_lv,min 

_lv,sec.lv]=Dataanalyze( 'test.lvm' ) ; 

4 

data_lv(:,1)=data 

_lv(: , 1) +str2num (sec.lv) +str2num (min.lv)*60+. .. 


str2num (h.lv)* 

3600; 

5 

%data_omni=csvread('example.txt'); 

6 

%data_omni (:,10) = 

12; 

7 

%data_omni ( :,14) = 

data.omni(:,13)/1000+data.omni(:,12)+data_omni... 


(:,11)*60+data_omni(:,10)*3600; 

8 

subplot (4,1,1); 


9 

plot (data.lv(:, 1) 

, data.lv ( : , 2) ) ; 

10 

subplot (4,1,2); 


11 

plot (data.lv(:,1) 

, data.lv ( : , 3) ) ; 

12 

subplot (4,1,3); 


13 

plot (data.lv(: , 1) 

, data.lv ( : , 4) ) ; 

14 

subplot (4,1,4); 


15 

plot (data.lv(: , 1) 

, data.lv ( : , 5) ) ; 

16 

%subplot(5,1,5); 


17 

%plot(data.lv(:,1 

) , data.lv ( :, 3) ) ; 

18 

%hold on 


19 

%plot(data.omni(: 

, 1) , data.lv(:,14)); 

B.1.2 Costume correlation function 1 


1 function [corr.out]=costume_corr(input_data) 

2 %input is [m,n]; m: length of each arr, n=numbers of arr. 

3 [-pn] =size (input.data); 

4 w=diag (ones(1,n)); 
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5 corr_out=zeros (n); 

6 for mm=1:n 

7 for mm2=mm:n 

8 corr_out(mm,mm2)=nancorrcoef(input-data(:,mm),... 

input-data(:,mm2)); 

9 end 

10 end 

u corr_out=corr_out + flipud(rot 90(corr_out))-w; 

12 end 


B.1.3 Costume correlation function 2 


1 clc 

2 close all 

3 clear all 

4 load('matlab.mat') 

5 load (' matlab_data.mat' ) 

6 load('matlab_numbers.mat') 

7 Total_finaldata=[Finaldata_nmtdayl,Finaldata_NMTday2,... 

FinaldataUNMday2]; 

8 Total_f inalNumb= [Finaldata_nmtdayl, Finaldata_NMTday2 , . . . 

FinaldataUNMday2]; 

9 list-load 

10 %% 


clc 

conver=[1 2 6 9 3 4 5 7 8 10 11 13 12 14 16 18 15 17 19 21 22 23 


20 6 11 3 10 13 22 2 
24 12 24 24 14 18]; 
ii=[1:125]; 


20 1 23 9 19 24 7 16 24 21 24 5 24 24 
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41 
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28 

31 

35 

27 

28 

31] ; 


15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


for d=l: length (yy) 
i=ii(d); 

nn=Total_finaldata(i).Name ; 
for qq=2: length (BoneLettering) 
e=BoneLettering(qq); 
w=e {1}; 

if (strncmpi(w,nn, min ([ length (w) length (nn)]))) 
test=l; 

match_numb=qq; 

end 

end 
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26 


y=conver(yy(d) ) ; 

27 if (Total_finaldata (i) .Name ( 1 ) =='A' ) 

28 Total_finaldata (i) .HoleType=l; 

29 elseif (Total_f inaldata (i).Name (1 )==' B ' ) 

30 Total.finaldata (i) .HoleType= 2 ; 

31 elseif (Total_f inaldata (i).Name (1 )==' C ' ) 

32 Total.finaldata (i) .HoleType= 3 ; 

33 elseif (Total_f inaldata (i).Name ( 1 ) == ' D ' ) 

34 Total_finaldata (i) .HoleType= 4 ; 

35 elseif (Total_f inaldata (i).Name (1 )==' E ' ) 

36 Total.finaldata (i) .HoleType= 5 ; 

37 end 

38 Total_f inaldata (i) . ID_UNMname=BoneLettering (matchmumb) ; 

39 Total_f inaldata (i) . DistancetoMark=DistancetoMark (matchmumb) ; 

40 Total_f inaldata (i) . FinalMa jorAxis=FinalMa jorAxis (matchmumb) ; 

41 Total_f inaldata (i) . FinalMinorAxis=FinalMinorAxis (matchmumb) ; 

42 Total_finaldata(i).InitialMajorAxisBoneplane= ... 

InitialMajorAxisBoneplane(matchmumb); 

43 Total_finaldata(i).InitialMinorAxisPerpendiculartoBone=... 

InitialMinorAxisPerpendiculartoBone(matchmumb); 

44 Total_f inaldata (i) . ID_number=form. ParticipantNumber_ 01 d (y) ; 

45 Total_f inaldata (i) . Years_in_residencty=form. . . . 

Years_in_Residency(y) ; 

46 Total_finaldata (i) . Sex=form. sex (y) ; 

47 Total_finaldata(i).Height=form.Height(y); 

48 Total_f inaldata (i) . Sex=form. sex (y) ; 

49 Total_f inaldata (i) . How_many_holes = form. How_many_holes (y) ; 

so Tot a l_f inaldata (i) . days_since_your_last_bone_drilling=form. . . 

days_since_your_last_bone_drilling (y) ; 

51 Total_f inaldata (i) . residency_program_include=form. . . . 

residency_program_include (y); 

52 Total_finaldata (i) . . . . 

program_have_dedicated_surgical_skills_training=form. . . . 
program_have_dedicated_surgical_skills_training (y) ; 

53 Total_f inaldata (i) . PVC=form. PVC (y) ; 

54 Total_f inaldata (i) . have_extensive_On_surgical_experience=form 

. have_extensive_On_surgical_experience (y) ; 

55 Total_f inaldata (i) . rate_your_drill_overpenetration=form. . . . 

rate_your_drill_overpenetration (y) ; 

56 %Total_finaldata (i) . rate_your_drill_overpenetration=form. . . . 

rate_your_drill_overpenetration (y) ; 

57 if (i< 41 ) 

58 Total_f inaldata (i) . Time_lv_ 2 ndcliped=Total_f inaldata (i) . . 

C_drilling_dataLV(:, 1 ); 

59 Total_f inaldata (i) . Force_lv_ 2 ndcliped=Total_f inaldata (i) . 

C_drilling_dataLV(:, 2 ); 

60 Total_f inaldata (i) . RPM _1 v_ 2 ndcliped=Total_f inaldata (i) . . . 

C_drilling_dataLV(:, 3 ); 

6 1 Total_finaldata (i) .bar_lv_ 2 ndcliped=Total_finaldata (i) . . . 

C_drilling_dataLV(:, 4 ); 
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62 


Total_finaldata(i).drill_lv_2ndcliped=Total_finaldata(i) 
C_drilling_dataLV(:,5); 

63 elseif (i<112) 

64 Total_f inaldata (i) . Time_lv_2ndcliped=Total_f inaldata (i) . 

C_drilling_dataLV(:,1); 

65 Total_f inaldata (i) . Force_lv_2ndcliped=Total_f inaldata (i) 

C_drilling_dataLV(:,4); 

66 Total_f inaldata (i) . RPM_lv_2ndcliped=Total_f inaldata (i) . . 

C_drilling_dataLV(:,5); 

67 Total_finaldata (i) .bar_lv_2ndcliped=Total_finaldata (i) . . 

C_drilling_dataLV(:,2); 

68 Total_f inaldata (i) . drill_lv_2ndcliped=Total_finaldata (i) 

C_drilling_dataLV(:,3); 

69 else 

70 Total_f inaldata (i) . Time_lv_2ndcliped=Total_f inaldata (i) . 

C_drilling_dataLV(:,1); 

71 Total_f inaldata (i) . Force_lv_2ndcliped=Total_f inaldata (i) 

C_drilling_dataLV(:,4); 

72 Tot a l_f inaldata (i) . RPM_1 v_2 n dcliped=Tot a l_f inaldata (i) . . 

C_drilling_dataLV(:,2); 

73 Total_finaldata (i) .bar_lv_2ndcliped=Total_finaldata (i) . . 

C_drilling_dataLV(:,5); 

74 Total_f inaldata (i) . drill_lv_2ndcliped=Total_f inaldata (i) 

C_drilling_dataLV(:,3); 

75 end 

76 [ res ] =rpm_f inal (Tot a l_f inaldata (i) . RPM_lv_2ndcl iped, . . . 

Total_finaldata (i) .Time_lv_2ndcliped); 

77 Total_f inaldata (i) . RPM_ave_2ndcliped=mean (res (1, : ) ) ; 

78 [w,-i] =size (Tot a l_f inaldata (i) . C_Overpen_fixed_dataMT) ; 

79 Total_finaldata (i) ,maxpen=sqrt ( (Total_finaldata (i) . . . . 

C_Overpen_fixed_dataMT(1,27)-Total_finaldata(i). ... 
C_Overpen_fixed_dataMT(w,27))"2+(Total_finaldata(i).... 
C_Overpen_fixed_dataMT(1,28)-Total_finaldata(i).... 
C_Overpen_fixed_dataMT(w,28))"2+(Total_finaldata(i).... 
C_Overpen_fixed_dataMT(1,29)-Total_finaldata(i)... . 
C_Overpen_fixed_dataMT(w,29))"2); 

so %Total_f inaldata (i) . avg_RPM=mean (Total_finaldata (i) . C_RPM. . . 

(:,2)) ; 

si if (Total_f inaldata (i) . RPM_ave_2ndcliped==0) 

82 Total_f inaldata (i) . RPM_ave_2ndcliped=NaN; 

83 end 

84 Total_f inaldata (i) . avg_force=mean (Total_f inaldata (i) .... 

Force_lv_2ndcliped) ; 

85 Total_finaldata (i) .max_f orce=max (Total_finaldata (i) .... 

Force_lv_2ndcliped) ; 

86 Total_finaldata (i) . std_f orce=std (Total_finaldata (i) .... 

Force_lv_2ndcliped) ; 

87 Total_f inaldata (i) . std_bar_vib=std (Total_finaldata (i) .... 

bar_lv_2ndcliped) ; 
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90 

91 

92 

93 


94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
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114 

115 

116 

117 

118 

119 

120 
121 


Total_f inaldata (i) . std_drill_vib=std (Total.f inaldata (i) .... 
drill_lv_2ndcliped) ; 

Total_finaldata (i) .tog_drill=nanmean([nanstd(Total_finaldata( 
i) . C_dr illing.f ixecLdat aMT (: , 21 ) ) , nanstd (Total_f inaldata ( . 
i) . C_drilling_f ixed_dataMT ( :, 2 8) ) , nanstd (T ot a l_f inaldata ( . 
i) . C_drilling_f ixed_dataMT (:, 29) ) ] ) ; 

Total.finaldata(i).mean_roll_drill=nanmean(Total_finaldata(i) 

. C_drilling_f ixed_dataMT ( : , 30) ) ; 

Total.finaldata (i) .mean.pitch_drill=nanmean(Total.finaldata(i 
) . C_drilling.fixed_dataMT ( :, 31) ) ; 

Total.finaldata(i).mean_yaw_drill=nanmean(Total.finaldata(i). 

C.drilling.fixed.dataMT ( :, 32)); 

%plot3(Total.finaldata(i) .C.drilling.fixed.dataMT (:,27 ), ... 
Total.f inaldata (i) . C.drilling.f ixed.dat aMT ( :, 28) , . . . 
Total.f inaldata (i) . C.drilling.f ixed.dat aMT ( :, 29) , 1 . 1 ) ; . . . 
axis equal; grid on 
%im=getframe; 

%imwrite(im.cdata,strcat(Total.finaldata(i).Name,’.JPG’)); 

end 

for i=112:125 

Total.finaldata(i).std_drill_vib=Total_finaldata(i).... 
std_drill_vib/10; 

end 

o_ o. 
o o 

Total.finaldata(51).Years_in_residencty=10; 

Total.finaldata(53) .Year s_in_residencty=l0; 

Total.finaldata (54) .Years_in_residencty=10; 

Total.finaldata (67) .Years_in_residencty=10; 

Total.finaldata(69).Years_in_residencty=l0; 

Total.finaldata(70).Years_in_residencty=10; 

Total.finaldata(81) .Year s_in_residencty=l0; 

Total.finaldata(83).Years_in_residencty=l0; 

Total.finaldata (84) .Years_in_residencty=l0; 

Total.finaldata(96).Years_in_residencty=l0; 

Total.finaldata(98) .Year s_in_residencty=l0; 

Total.finaldata(99) .Year s_in_residencty=l0; 

Total.finaldata(109) .Years _in_residencty=10; 

Total.finaldata(110) .Year s.in.residencty=10; 

Total.finaldata(114) .Year s.in.residencty=10; 

Total.finaldata(115) .Year s.in.residencty=10; 

Total.finaldata(117) .Year s_in_residencty=10; 

o, o. 
o o 

finalp=costume_corr([Total.finaldata.HoleType;Total.finaldata .... 
DistancetoMark;Total.finaldata.FinalMajorAxis;Total.finaldata. 
.FinalMinorAxis;Total.finaldata.ID.number;Total.finaldata.... 
Years.in.residencty;Total.finaldata.Sex;Total.finaldata.... 
Height;Total.finaldata.How_many_holes; Total.finaldata. ... 
days.since.your.last.bone.drilling; Total.finaldata. ... 
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residency_program_include; Total_finaldata. . . . 
program_have_dedicated_surgical_skills_training; . . . 

Total_finaldata.PVC; Total_finaldata. ... 

have_extensive_On_surgical_experience; Total_f inaldata . . . . 
RPM_ave_2ndcliped;Total_finaldata.maxpen;Total_finaldata.... 
avg_force;Total_finaldata.max_force;Total_finaldata.std_force 
;Total_finaldata.std_bar_vib;Total_finaldata.std_drill_vib;.. 
Total_finaldata.tog_drill;Total_finaldata.mean_roll_drill;... 
Total_finaldata.mean_pitch_drill;Total_finaldata.... 
mean_yaw_drill]'); 

122 % l-Total_f inaldata . HoleType; 

123 % 2-Total_f inaldata . DistancetoMark; 

124 % 3-Total_f inaldata . FinalMajorAxis; 

125 % 4-Total_f inaldata .FinalMinorAxis; 

126 % 5-Total_f inaldata . ID.number; 

127 % 6-Total_f inaldata. Years_in_residencty; 

128 % 7-Total_f inaldata . Sex; 

129 % 8-Total_f inaldata . Height ; 

130 % 9-Total_f inaldata . How_many_holes ; 

131 % 10-Total_f inaldata . days_since_your_last_bone_drilling; 

132 % ll-Total_f inaldata . residency_program_include; 

133 % 12-Total_f inaldata . . . . 

program_have_dedi cat ed_surgi cal _skills_t raining; 

134 % 13-Total_f inaldata . PVC; 

135 % 14-Total_f inaldata . have_extensive_On_surgical_experience; 

136 % 15-Total_f inaldata . RPM_ave_2ndcliped; 

137 % 16-Total_f inaldata . maxpen; 

us % 17-Total_finaldata.avg_force; 

139 % 18-Total_f inaldata . max_f orce ; 

140 % 19-Total_f inaldata . std_f orce ; 

141 % 20-Total_f inaldata . std_bar_vib; 

142 % 21-Total_f inaldata . std_drill_vib; 

143 % 22-Total_f inaldata . tog_drill; 

144 % 23-Total_f inaldata . mean_roll_drill; 

145 % 24-Total_f inaldata . mean_pitch_drill ; 

146 % 25-Total_f inaldata . mean_yaw_drill ; 

147 imagesc ((f inalp) ) %plot the 

148 %% 

149 figure 

150 pre_Total_finaldata=Total_finaldata (1:40) ; 

151 figure; pre_f inalp=costume_corr ( [ pre_Tot a l_f inaldata . HoleType; . . . 

pre_Total_f inaldata . DistancetoMark; pre_Total_f inaldata . . . . 
FinalMajorAxis;pre_Total_finaldata.FinalMinorAxis; . . . 
pre_Total_finaldata.ID_number;pre_Total_finaldata . . . . 
Years_in_residencty; pre_Total_f inaldata . Sex; . . . 
pre_Total_f inaldata . Height; pre_Total_f inaldata . How_many_holes 
; pre_Total_f inaldata . days_since_your_last_bone_drilling; ... 
pre_Total_f inaldata . residency_program_include; . . . 
pre_Total_finaldata . . . . 

program_have_dedicated_surgical_skills_training; . . . 
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pre_Total_finaldata.PVC; pre_Total_finaldata.... 
have_extensive_On_surgical_experience; pre_Total_f inaldata . . . . 
RPM_ave_2ndcliped; pre_Total_f inaldata .maxpen; . . . 
pre_Total_f inaldata . avg.force; pre_Total_f inaldata . max_f orce; . 
pre_Total_f inaldata . std_force; pre_Total_f inaldata . std_bar_vib 
; pre_Total_f inaldata . std_drill_vib; pre_Total_f inaldata . . . . 
tog_drill; pre_Total_f inaldata . mean.rol l_dril 1; . , , 
pre_Total_f inaldata . mean_pitch_drill; pre_Total_f inaldata . . . . 
mean_yaw_drill]'); 


152 imagesc ( (pre_f inalp) ) 

153 %title ( ' PRE ' ) 

154 post_Total_finaldata=Total_finaldata ( 41:12 5 ) ; 

155 figure; Post_f inalp=costume_corr ( [post_Total_f inaldata . HoleType; . 

post_Total_finaldata.DistancetoMark;post_Total_finaldata. . . . 
FinalMajorAxis;post_Total_f inaldata.FinalMinorAxis; . . . 
post_Total_finaldata.ID.number;post_Total_finaldata. ... 
Years_in_residencty; post_Total_f inaldata . Sex; . . . 
post_Total_f inaldata.Height;post_Total_finaldata . . . . 
How_many_holes; post_Total_finaldata . . . . 

days_since_your_last_bone_drilling; post_Total_finaldata . . . . 
residency_program_include; post_Total_finaldata . . . . 
program_have_dedicated_surgical_skills_training; . . . 
post_Total_finaldata .PVC; post_Total_finaldata . . . . 
have_extensive_On_surgical_experience; post_Total_finaldata . . . 
RPM_ave_ 2 ndcliped; post_Total_f inaldata .maxpen; . . . 
post_Total_f inaldata . avg.force; post_Total_f inaldata .max.force 
;post_Total_f inaldata.std.force;post_Total_finaldata.... 
std_bar_vib;post_Total_finaldata.std_drill_vib;... 
post_Total_f inaldata . tog_drill; post_Total_f inaldata . ... 
mean_roll_dr ill; post_Total_f inaldata . mean_pitch_dri 11 ; . . . 
post_Total_f inaldata . mean_yaw_drill ] ' ) ; 

156 imagesc ( (Post_f inalp) ) 

157 title ( ' POST ' ) 

158 diff_finalp=Post_finalp-pre_finalp; 

159 figure 

wo imagesc ((diff_finalp)) 
i 6 i title (' dif f_f inalp ' ) 


162 


165 


163 


164 


figure 
hold on 
for i=l:125 


172 


166 


167 


168 


169 


170 


171 


if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
DistancetoMark,' c*' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
DistancetoMark, '*m') 
case 3 
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173 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
DistancetoMark, 'b*' ) 

174 case 4 

175 plot (Total_finaldata (i) . ID_number, Total_finaldata (i) 

DistancetoMark, 'r*' ) 

176 case 5 

177 plot (Total_finaldata (i) . ID.number, Total_finaldata (i) 

DistancetoMark, 'g*' ) 
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184 

185 

186 

187 

188 

189 

190 

191 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
DistancetoMark,' .c' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
DistancetoMark,' .m' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
DistancetoMark,' .b' ) 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
DistancetoMark,' .r' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
DistancetoMark,' .g' ) 
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195 
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199 

200 
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204 

205 

206 

207 

208 

209 

210 

211 


end 

end 

end 

title ( 'Distanc2marker' ) 
im=getf rame ; 

imwrite(im.cdata,strcat( 'Distanc2marker' , ' .JPG' )) ; 

figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis, 'c*' ) 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMajorAxis, 'm*' ) 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis, 'b*' ) 
case 4 
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212 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis, 'r*' ) 

213 case 5 

214 plot (Total_finaldata (i) . ID_number, Total_finaldata (i) 

FinalMajorAxis, 'g*' ) 


215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis,' .c' ) 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMajorAxis,' .m' ) 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis,' .b' ) 
case 4 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMajorAxis,' .r' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMajorAxis,' .g' ) 


229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 


end 

end 

end 

title ( 'FinalMajorAxis' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'FinalMajorAxis' ,' .JPG' )); 

figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis, 'c*' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis, 'm*' ) 
case 3 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
FinalMinorAxis, 'b*' ) 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis, 'r*' ) 
case 5 
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251 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis,' g*' ) 


252 

253 

254 

255 

256 

257 

258 

259 

260 
261 

262 

263 

264 

265 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMinorAxis,' .c' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis,' .m' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMinorAxis,' .b' ) 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
FinalMinorAxis,' .r' ) 
case 5 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
FinalMinorAxis,' .g' ) 


266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 
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278 
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280 
281 
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285 
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287 


end 

end 

end 

title ( 'FinalMinorAxis' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'FinalMinorAxis' ,' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
RPM_ave_2ndcliped, ' c* ' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
RPM_ave_2ndcliped, 'm* ' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
RPM_ave_2ndcliped, 'b*') 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
RPM_ave_2ndcliped, 'r*' ) 
case 5 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
RPM_ave_2ndcliped, 'g*') 


288 

289 end 
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290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 


else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
RPM_ave_2ndcliped,' .c' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
RPM_ave_2ndcliped, ' .m' ) 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
RPM_ave_2ndcliped, ' .b' ) 
case 4 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
RPM_ave_2ndcliped, '.r' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
RPM_ave_2ndcliped, ' .g' ) 
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305 
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307 
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309 
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315 

316 

317 

318 

319 

320 

321 
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323 


end 

end 

end 

title ( ' RPM_ave_2ndcliped' ) 
im=getf rame ; 

imwrite(im.cdata,strcat( 'RPM_ave_2ndcliped' ,' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
maxpen, 'c*' ) 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
maxpen, 'm*' ) 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
maxpen,' b*' ) 
case 4 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
maxpen,' r*' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
maxpen, 'g*' ) 


324 

325 

326 

327 

328 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 
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337 
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339 
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342 
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359 
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361 

362 

363 

364 

365 

366 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
maxpen, '.c' ) 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
maxpen, '.m' ) 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
maxpen, '.b' ) 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
maxpen, '.r' ) 
case 5 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
maxpen, '.g' ) 


end 

end 

end 

title ( 'maxpen' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'maxpen' ,' . JPG' ) ) ; 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
avg_force, 'c*' ) 
case 2 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
avg.force, 'm*') 
case 3 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
avg_force, 'b*' ) 
case 4 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
avg_force, ' r* ' ) 
case 5 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
avg_force, ' g* ' ) 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
avg_force,' .c' ) 
case 2 
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367 


368 

369 

370 

371 

372 

373 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
avg_force, ' .m' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
avg_force, ' .b' ) 
case 4 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
avg_force, ' .r' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
avg_force, ' .g' ) 
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388 

389 

390 

391 

392 

393 
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395 


end 

end 


end 

title ( 'avg_force' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'avg_force' , '.JPG' )); 

figure 

hold on 

for i = l : 125 


if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
std_bar_vib, ' c* ' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
std_bar_vib, 'm* ' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
std_bar_vib, 'b*') 
case 4 

plot (Total_finaldata (i) . IDmumber, Total_finaldata (i) 
std_bar_vib, ' r* ' ) 
case 5 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
std_bar_vib, ' g* ' ) 


396 

397 

398 

399 

400 

401 

402 

403 

404 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
std_bar_vib, ' . c ' ) 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
std_bar_vib, ' .m' ) 
case 3 
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405 


plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
std_bar_vib, ' .b' ) 

406 case 4 

407 plot (Total_finaldata (i) . ID_number, Total_finaldata (i) 

std_bar_vib, ' . r ' ) 

408 case 5 

409 plot (Total_finaldata (i) . ID.number, Total_finaldata (i) 

std_bar_vib, ' .g' ) 
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414 
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416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 


end 

end 

end 

title ( ' std_bar_vib ' ) 
im=getf rame ; 

imwrite(im.cdata,strcat( 'st d_bar_vib' , ' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
mean_roll_drill, 'c*') 
case 2 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
mean_roll_drill, 'm*') 
case 3 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
mean_roll_drill, 'b*') 
case 4 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
mean_roll_drill, 'r*') 
case 5 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
mean_roll_drill, 'g*') 


432 
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435 

436 

437 

438 

439 

440 

441 

442 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
mean_roll_drill,' .c' ) 
case 2 

plot (Total_finaldata(i).ID.number,Total_finaldata(i) 
mean_roll_drill,' ,m' ) 
case 3 

plot (Total_finaldata(i) .ID_number,Total_finaldata (i) 
mean_roll_drill,' ,b' ) 
case 4 
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443 plot (Total_finaldata (i) . ID.number , Total_finaldata (i) . 

mean_roil-drill, ' . r ' ) 

444 case 5 

445 plot (Total_finaldata (i) . ID.number, Total_finaldata (i) . 

mean_roll_drill , ' .g' ) 
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end 

end 

end 

title ( ’ mean_roll_drill 1 ) 
im=getframe; 

imwrite(im.cdata,strcat( 1 mean.roll_dri11 ' , ? .JPG’ )); 


2 - 2 - 
o o 

figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

DistancetoMark, 'c*' ) 
case 2 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

DistancetoMark, '*m' ) 
case 3 

plot (Total-finaldata(i).HoleType,Total_finaldata(i). . 

DistancetoMark, 'b*' ) 
case 4 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

DistancetoMark, 'r*' ) 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 
DistancetoMark, 'g*' ) 


471 
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474 

475 

476 

477 

478 

479 
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481 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
.DistancetoMark,' .c' ) 
case 2 

plot (Total-finaldata(i).HoleTypet5,Total_finaldata(i) 
.DistancetoMark,' ,m' ) 
case 3 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
•DistancetoMark,' .b' ) 
case 4 
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482 


plot (Total.finaldata(i).HoleType+5,Total-finaldata(i) 
.DistancetoMark,' .r' ) 

483 case 5 

484 plot (Total_f inaldata (i) . HoleType+5, Total_f inaldata (i) 

.DistancetoMark,' .g' ) 
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end 

end 

end 

title ( 'Distanc2marker' ) 
im=getf rame ; 

imwrite(im.cdata,strcat( 'Distanc2marker+5' , ' .JPG' )); 

figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total-finaldata(i).HoleType,Total_finaldata(1).. 

FinalMajorAxis, 'c*' ) 
case 2 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

FinalMajorAxis, 'm*' ) 
case 3 

plot (Total-finaldata(i).HoleType,Total_finaldata(1).. 

FinalMajorAxis, 'b*' ) 
case 4 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

FinalMajorAxis, 'r*' ) 
case 5 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 
FinalMajorAxis, 'g*' ) 
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520 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMajorAxis,'. c' ) 
case 2 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMajorAxis,' .m' ) 
case 3 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMajorAxis,' .b' ) 
case 4 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMajorAxis,' .r' ) 
case 5 
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559 


plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMajorAxis, '.g' ) 


end 

end 

end 

title ( 'FinalMajorAxis' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'FinalMajorAxis+5' ,' .JPG' )); 

figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total.finaldata(i).HoleType,Total_finaldata(1).. 

FinalMinorAxis, 'c*' ) 
case 2 

plot (Total-finaldata(i).HoleType,Total_finaldata(1).. 

FinalMinorAxis, 'm*' ) 
case 3 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

FinalMinorAxis, 'b*' ) 
case 4 

plot (Total-finaldata(i).HoleType,Total_finaldata(1).. 

FinalMinorAxis, 'r*' ) 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 
FinalMinorAxis, 'g*' ) 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMinorAxis,'. c' ) 
case 2 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMinorAxis,' .m' ) 
case 3 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMinorAxis,' .b’ ) 
case 4 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMinorAxis,' .r' ) 
case 5 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.FinalMinorAxis,' .g' ) 


68 




560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 


end 

end 

end 

title ( 'FinalMinorAxis' ) 
im=get frame; 

imwrite(im.cdata,strcat ( ' FinalMinorAxis+5 ' ,' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

RPM_ave_2ndcliped, 'c*' ) 
case 2 

plot (Total.finaldata(i).HoleType,Total_finaldata(i).. 

RPM_ave_2ndcliped, 'm* ' ) 
case 3 

plot (Total-finaldata(i).HoleType,Total_finaldata(i). . 

RPM_ave_2ndcliped, 'b* ' ) 
case 4 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

RPM_ave_2ndcliped, 'r*' ) 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 
RPM_ave_2ndcliped, ' g* ' ) 
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end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. RPM_ave_2ndcliped, ' . c ' ) 
case 2 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
. RPM_ave_2ndcliped, ' .m') 
case 3 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. RPM_ave_2ndcliped, ' . b ' ) 
case 4 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. RPM_ave_2ndcliped, ' . r ' ) 
case 5 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
. RPM-ave_2ndcliped, ' . g' ) 


595 

596 end 

597 end 

598 end 

599 title (' RPM-ave_2ndcliped ' ) 
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620 
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im=getf rame ; 

imwrite(im.cdata,strcat( 'RPM_ave_2ndcliped+5' ,' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

maxpen, 'c*' ) 
case 2 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

maxpen, 'm*' ) 
case 3 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

maxpen,' b*' ) 
case 4 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

maxpen,' r*' ) 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i). . 
maxpen, 'g*' ) 

end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
.maxpen, '.c' ) 
case 2 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
.maxpen, '.m' ) 
case 3 

plot (Total-finaldata(i).HoleTypet5,Total_finaldata(i) 
.maxpen, '.b' ) 
case 4 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
.maxpen, '.r' ) 
case 5 

plot (Total_finaldata(i).HoleTypet5,Total_finaldata(i) 
.maxpen, '.g' ) 

end 

end 

end 

title ( 'maxpen' ) 
im=getf rame ; 

imwrite(im.cdata,strcat( 'maxpen+5' ,' .JPG' )); 

figure 

hold on 
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for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total.finaldata(i).HoleType,Total_finaldata(i).. 

avg.force, ' c* ' ) 
case 2 

plot (Total-finaldata(i).HoleType,Total_finaldata(i). . 

avg.force, 'm*') 
case 3 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

avg.force, 'b*' ) 
case 4 

plot (Total_finaldata(i).HoleType,Total_finaldata(i).. 

avg_force, 'r*' ) 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 
avg_force, 'g*' ) 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
.avg_force,' .c' ) 
case 2 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.avg_force,' .m' ) 
case 3 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.avg_force,' .b' ) 
case 4 

plot (Total-finaldata(i).HoleType+5,Total_finaldata(i) 
.avg-force, '.r' ) 
case 5 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.avg-force,' .g' ) 


667 

668 end 

669 end 

670 end 

671 title (' avg-force ' ) 

672 im=getf rame ; 

673 imwrite (im. cdata, strcat ( ' avg_force + 5 ' , ' . JPG ' ) ) ; 

674 figure 

675 hold on 

676 for i=l: 125 

677 if (i<40) 

678 switch Total_f inaldata (i) . HoleType 

679 case 1 
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680 


681 

682 

683 

684 

685 

686 

687 

688 


plot (Total_finaldata(i) 
stcLbar_vib, 'c*' ) 
case 2 

plot (Total_finaldata(i) 
std_bar_vib, 'm* ' ) 
case 3 

plot (Total_finaldata(i) 
std_bar_vib, 'b* ' ) 
case 4 

plot (Total_finaldata(i) 
std_bar_vib, ' r* ' ) 
case 5 

plot (Total_finaldata(i) 

std_bar_vib, ' g* ' ) 


.HoleType,Total_finaldata(i).. 


.HoleType,Total_finaldata(i).. 

.HoleType,Total_finaldata(i).. 

.HoleType,Total_finaldata(i).. 


.HoleType,Total_finaldata(i).. 
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end 


else 


switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. std_bar_vib, ' . c ' ) 
case 2 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.std_bar_vib, ' .m' ) 
case 3 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. std_bar_vib, ' . b ' ) 
case 4 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. std_bar_vib, ' . r ' ) 
case 5 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
. std_bar_vib, ' . g' ) 
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end 

end 

end 

title ( ' std_bar_vib ' ) 
im=get frame; 

imwrite(im.cdata,strcat( 'st d_bar_vib + 5' , ' .JPG' )); 
figure 
hold on 
for i=l:125 
if (i<4 0) 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

mean_roll_drill, 'c*') 
case 2 
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plot (Total.finaldata(i).HoleType,Total_finaldata(i).. 

mean_roll_drill, 'm*') 
case 3 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 

mean_roll_drill, 'b*') 
case 4 

plot (Total-finaldata(i).HoleType,Total_finaldata(1).. 

mean_roll_drill, 'r*') 
case 5 

plot (Total-finaldata(i).HoleType,Total_finaldata(i).. 
mean_roll_drill, 'g*') 


end 

else 

switch Total_finaldata(i).HoleType 
case 1 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.mean_roll_dri11, '.c' ) 
case 2 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.mean_roll_drill, '.m' ) 
case 3 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.mean_roll_drill, '.b' ) 
case 4 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.mean.roll_dri11, '.r' ) 
case 5 

plot (Total_finaldata(i).HoleType+5,Total_finaldata(i) 
.mean_roll_dri11, '.g' ) 


739 

740 end 

741 end 

742 end 

743 title ( 'mean-roll-drill ' ) 

744 im=getf rame ; 

745 imwrite (im. cdata, strcat ( 'mean_roll_drill + 5 ' , ' . JPG' ) ) ; 

746 %% 

747 pullout= [1:51; 2447 1809 2122 2447 2447 2876 2695 2780 2366 3075 

2409 2765 2568 2945 2932 nan 2423 2887 3398 2376 2371 2839 .. 

2715 2428 2145 2075 3089 3204 2626 2851 2525 2362 2787 2714 . 

2497 3096 2209 3185 2613 2985 3032 2825 1937 3528 2623 3409 . 

2933 2439 3341 3343 2905]'; 

748 plot (pullout ( : , 1) , pullout (: , 2) , ' . ' ) 

749 ii=l; 

750 for i=l: 125 

751 if (Total_finaldata (i) .HoleType==3) 

752 PulloutTotal_finaldata (ii) =Total_finaldata (i) ; 

753 PulloutTotal_finaldata (ii) .HoleType=3; 

754 ii=ii + l; 
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755 

end 

756 

end 

757 

2- 2- 
o o 

758 

2- 2- 
o o 

759 

figure 

760 

finalp— 


DistancetoMark;Total.finaldata.FinalMajorAxis;Total.finaldata 
.FinalMinorAxis;Total.finaldata.ID.number;Total.finaldata.... 
Years.in.residencty;Total.finaldata.Sex;Total.finaldata.... 
Height;Total.finaldata.How_many_holes; Total.finaldata. ... 
days.since.your.last.bone.drilling; Total.finaldata.. . . 
residency.program.include; Total.finaldata. . . . 
program.have.dedicated.surgical.skills.training; ... 

Total.finaldata.PVC; Total.finaldata. ... 

have.extensive.On.surgical.experience; Total.f inaldata . . . . 
RPM_ave_2ndcliped;Total.finaldata.maxpen;Total.finaldata.... 
avg.force;Total.finaldata.max.force;Total.finaldata.std.force 
;Total.finaldata.std.bar.vib;Total.finaldata.std.drill.vib;.. 
Total.finaldata.tog.drill;Total.finaldata.mean.roll.drill;... 
Total.finaldata.mean.pitch.dri11;Total.finaldata.... 
mean.yaw.drill]'); 

761 imagesc ((f inalp) ) %plot the 

762 % l-Total_f inaldata . HoleType; 

763 % 2-Total.f inaldata . DistancetoMark; 

764 % 3-Total_f inaldata .FinalMa jorAxis; 

765 % 4-Total_f inaldata . FinalMinorAxis; 

766 % 5-Total_finaldata . ID.number ; 

767 % 6-Total_f inaldata . Years.in.residencty; 

768 % 7-Total.f inaldata . Sex; 

769 % 8-Total.f inaldata . Height ; 

770 % 9-Total_f inaldata . How.many.holes ; 

771 % 10-Total_f inaldata . days.since.your.last.bone.drilling; 

772 % ll-Total_f inaldata . residency.program.include; 

773 % 12-Total.f inaldata . . . . 

program.have.dedi cat ed.surgi cal .ski lls.t raining; 

774 % 13-Total_f inaldata . PVC; 

775 % 14-Total_f inaldata . have.extensive.On.surgical.experience; 

776 % 15-Total_f inaldata . RPM_ave_2ndcliped; 

777 % 1 6-Total_f inaldata . maxpen; 

778 % 17-Total_f inaldata . avg.force; 

779 % 18-Total_f inaldata . max.f orce ; 

780 % 19-Total_f inaldata . std.f orce ; 

781 % 20-Total_f inaldata . std.bar.vib; 

782 % 21-Total_f inaldata . std.drill.vib; 

783 % 22-Total_f inaldata . tog.drill ; 

784 % 23-Total_f inaldata . mean.roll.drill ; 

785 % 24-Total_f inaldata . mean.pitch.drill ; 

786 % 25-Total_f inaldata . mean.yaw.drill ; 

787 

788 %% 
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789 figure 

790 for i=l: length (PulloutTotal_finaldata) 

791 PulloutTotal_finaldata(i) . ID_number 

792 PulloutTotal_finaldata (i) .Pullout=pullout ( . . . 

PulloutTotal_finaldata(i).ID.number,2); 

793 end 

794 %Pul lout Tot a l_f inaldata . HoleType ; 

795 f inalp=costume_corr ( [PulloutTotal_f inaldata . DistancetoMark; . . . 

PulloutTotal_finaldata.FinalMajorAxis;PulloutTotal_finaldata.... 
FinalMinorAxis;PulloutTotal_finaldata.ID_number;... 
PulloutTotal_finaldata.Years_in_residencty; . . . 

PulloutTotal_finaldata.Sex;PulloutTotal_finaldata.Height;... 
PulloutTotal_finaldata.How_many_holes; PulloutTotal_finaldata... 
.days_since_your_last_bone_drilling; PulloutTotal_finaldata. ... 
residency_program_include; PulloutTotal_finaldata.... 
program_have_dedicated_surgical_skills_training; . . . 
PulloutTotal_finaldata.PVC; PulloutTotal_finaldata. ... 
have_extensive_On_surgical_experience; PulloutTotal_finaldata . . . . 
RPM_ave_2ndcliped;PulloutTotal_finaldata.maxpen; . . . 
PulloutTotal_finaldata.avg_force; PulloutTotal_finaldata. ... 
max.force;PulloutTotal_finaldata.std.force;... 

PulloutTotal_finaldata.std_bar_vib;PulloutTotal_finaldata. ... 
std_drill_vib;PulloutTotal_finaldata.tog_drill; . . . 

Pu1loutTotal_finaldata.mean_roll_drill;PulloutTotal_finaldata... 
.mean_pitch_drill;PulloutTotal_finaldata.mean_yaw_dril1;... 
PulloutTotal.finaldata.Pullout] ') ; 

796 imagesc ((f inalp) ) 

797 im=get frame; 

798 imwrite(im.cdata,strcat( 'COV_with _pull_out ' , ' . JPG' ) ) ; 


B.1.4 Data fitter function 


1 function [a]=data_fit(lavbiefile, hapticfile) 

2 [dataLV_A,H,MIN, SEC]=Labview2var(lavbiefile); 

3 haptic_A=csvread (hapticfile); 

4 %% 

5 [dataLV_A(:,2)]=Low_filter.force(dataLV_A(:,2)); 

6 dataLV_A (: , 2)=5 601.7*dataLV_A (:,2)+6208.7; 

7 figure; plot (dataLV_A(: , 1),dataLV_A(:,2)); 

8 haptic_A (: , 3) =str2num (H); 

9 haptic_A (: , 7 ) =haptic_A (: , 3)*3600+haptic_A (: , 4)* 60+haptic_A(: , 5)+ . . . 

haptic_A (: , 6) *0.001; 

10 hold on; plot (hapt ic_A ( : , 7 ) , hapt ic_A ( : , 1) * 10, ' r ' ) ; plot (hapt ic_A . . . 

(:,7),haptic_A(:,2)*10000, 'g') 

u %% 

12 clip_A=[6.398e+04 1.464e+05]; 
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13 [OUTPUT] =syncer (dataLV_A, haptic_A, clip_A) ; 

14 figure 

is sf = fit([OUTPUT(:,1), OUTPUT(:,2)],OUTPUT(:,3),'poly23') 

16 plot(sf,[OUTPUT(:,1), OUTPUT (:, 2)], OUTPUT (:, 3) ) 

17 a=l; 
is end 


B.1.5 Row data adder function 


1 clc 

2 clear all 

3 % [data_lv,h_lv,min.lv,sec.lv]=Dataanalyze('test.lvm'); 

4 % data.lv(: , 1)=data.lv(: , 1)+str2num(sec.lv)+str2num(min.lv)*60+.. . 

str2num(h.lv)*3600; 

5 data_omni=csvread ( '2.txt' ); 

6 %data_omni(:,10)=12; 

7 data.omni(:,17)=data_omni(:,13)/1000+data_omni(:,12)fdata.omni ... 

(: , 11)*60+data_omni(:,10)*3600; 

8 for q=l: length (data.omni) 

9 [x,y,z]=trans_omni([data.omni(q,10),data.omni(q,11),data.omni... 

(q,12)], [data.omni(q,1),data.omni(q,2),data.omni(q, 3)], [ ... 
data.omni(q,7),data.omni(q,8),data.omni(q,9) ] ,[-174 -81 ... 

0] ) ; 

10 Location.drill(q, 1)=x; Location.drill(q, 2)=y; Location.drill(q. . . 

,3)=z; 

n end 

12 plot3 (Location.drill(: , 1),Location.drill(: , 2),Location.drill(:,3) . . . 

) 

13 grid on 

14 axis equal 


B.1.6 General analyzer function 


1 function [dataLV,H,MIN,SEC]=LTFO(filename) 

2 filelD = fopen (filename); 

3 fid=fileID; 

4 C = textscan (filelD, '%s %s'); 

5 Timeloc=ll; 

6 data_start=30; 

7 com=strfind(C{1,l}{Timeloc,1}, ',' ); 

8 colo=strfind(C{l,l}{Timeloc,1},':'); 

9 H= C{1,l}{Timeloc,1}((com+1):(colo(1,1)-1)); 
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10 MIN= C{1,l}{Timeloc,1}((colo(1,1)+1):(colo(1,2)-1)); 

11 SEC= C{l,l}{Timeloc,1}( (colo (1,2)+1) : (colo (1,2)+8)); 

12 fid=fopen (filename) ; 

13 f id2=f open ( ' newf ile . txt' , ' wt' ) ; 

14 id= 0 ; 

is a=f gets (fid) ; 

16 while (ischar (a) ) 

17 id=id+l; 

is if id<24 

19 a=fgets (fid) ; 

20 continue 

21 else 

22 fprintf (fid2, a) ; 

23 end 

24 a=fgets (fid) ; 

25 end 

26 dataLV=csvread ( ' newf ile . txt' ) ; 


B.1.7 Low pass filter function 1 


1 d = designfilt('lowpassfir'); 

2 y = filtfilt(d,data); 

3 yl = filter (d,data); 

4 subplot (2,1,1) 

5 plot ([y yl]) 

6 title('Filtered Waveforms') 

7 legend (' Zero-phase FilteringConventional Filtering') 

8 

9 subplot (2,1,2) 

10 plot (data) 

u title('Original Waveform') 


B.1.8 Main file 2 


1 clc 

2 clear all 

3 %List of data files to import: 

4 list=[ 'R54 ' ; 'R53 ';'R52 ' ; 'R43P' ; 'R43 ' ; 'R42P' ; 'R42 ' ; 'R4IP 

R41 ' ; 'R2P3' ; 'R2P2' ; 'R2P1' ; 'R23 ';'R22 ';'R21 ' ; 'R13P' ; 'R13 
; 'R12P' ;' R12 ' ; 'R11P' ;' Rll ';'MS43 ]; 

5 List=cellstr(list); 

6 Points=[96344 119417 120446 1; 
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7 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


32 

33 

34 

35 

36 


37 

38 

39 


40 


41 


96297 126726 127559 2; 

121617 138161 138161 3; 

66880 123987 124686 4; 

65301 150367 152001 5; 

84015 137039 139502 6; 

99682 162385 164541 7; 

112468 167029 168706 8; 

65276 122882 133961 9; 

57996 101332 104183 10; 

102031 113035 113617 11; 

65935 99735 100861 12; 

105174 123514 124200 13; 

57503 67468 70365 14; 

67010 88523 89537 15; 

111229 129450 130126 16; 

87093 126329 142184 17; 

75894 116583 117238 18; 

75426 105250 106975 19; 

86157 121960 122882 20; 

104253 151079 153191 21; 

83437 103411 105106 22]; 
for i = l: length (list); 

data(i).name = char (List (i,1)); 
data (i).test = Points (i,:); 

[data(i).NI,h_lv,min_lv,sec_lv]=Dataanalyze(strcat(char(List( 
i,1 )),' .lvm' )); 

data(i).NI(:,1)=data(i).NI(:,1) +str2num (sec_lv) +str2num (... 

min.lv) *60+str2num (h_lv)*3600; 
data(i). Omni=csvread (strcat(char(List(i,l)),' .txt' )); 
data(i).Omni(:,17)=data(i).Omni(:,13)/1000+data(i).Omni(:,12) 
+data(i).Omni(:,11)*60+data(i).Omni(:,10)*3600; 
for q=l: length (data(i).Omni) 

[x, y,z]=trans_omni([data(i) .Omni(q,10),data(i) .Omni(q,11) 

,data(i).Omni(q,12)],[data(i).Omni(q,1),data(i).Omni( . 
q,2),data(i).Omni(q,3)],[data(i).Omni(q,7),data(i).. .. 
Omni(q,8),data(i).Omni(q,9)],[-174 -81 0]); 
data (i) .Location(q,1)=x;data(i) .Location(q,2)=y;data(i) . . 
Location(q,3)=z; 

end 

data(i).Holelength=((data(i).Location(Points(i,1),1)-data(i). 
Location(Points(i,2),1))~2+(data(i).Location(Points(i,1).. 
,2)-data(i).Location(Points(i,2),2))~2 +(data(i).Location . 
(Points(i,1),3)-data(i).Location(Points(i,2),3))~2)".5; 
data(i) .Overpenetrationlength=((data(i) .Location (Points (i,2) . 
,1)-data(i).Location(Points(i,3),l))~2+(data(i).Location( . 
Points(i,2),2)-data(i).Location(Points(i,3),2))~2 +(data( . 
i) .Location(Points(i,2),3)-data(i) .Location(Points(i,3) . . . 
,3)) ~ 2 ) ~.5; 

for num=Points(i,1):Points(i,2) 
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42 


data(i).Penetrationspeed(num-Points(i,1)+1)=(data(i).Omni... 
(num,4)~2+data(i).Omni(num,5)~2+data(i).Omni(num,6 )... 

*2)~.5; 

43 end 

44 data(i).MaxPenetrationspeed=max (abs (data(i).Penetrationspeed... 

)); 

45 data(i).MinPenetrationspeed=min (abs (data(i).Penetrationspeed... 

)); 

46 end 


B.1.9 File loader (Customized Matlab generated code) 


1 function [BoneLettering,InitialMajorAxisBoneplane,... 

InitialMinorAxisPerpendiculartoBone,FinalMajorAxis, . . . 
FinalMinorAxis,DistancetoMark,NOTE] = importfile(workbookFile... 

,sheetName,startRow,endRow) 

2 %IMPORTFILE Import data from a spreadsheet 

3 % [BoneLettering,InitialMajorAxisBoneplane, ... 

InitialMinorAxisPerpendiculartoBone,FinalMajorAxis, . . . 
FinalMinorAxis,DistancetoMark,NOTE] 

4 % = IMPORTFILE(FILE) reads data from the first worksheet in the... 

Microsoft 

5 % Excel spreadsheet file named FILE and returns the data as ... 

column 

6 % vectors. 

7 % 

8 % [BoneLettering,InitialMajorAxisBoneplane, ... 

InitialMinorAxisPerpendiculartoBone,FinalMajorAxis, ... 
FinalMinorAxis,DistancetoMark,NOTE] 

9 % = IMPORTFILE(FILE,SHEET) reads from the specified worksheet. 

10 % 

u % [BoneLettering,InitialMajorAxisBoneplane,... 

InitialMinorAxisPerpendiculartoBone,FinalMa jorAxis,... 
FinalMinorAxis,DistancetoMark,NOTE] 

12 % = IMPORTFILE (FILE,SHEET,STARTROW, ENDROW) reads from the ... 

specified 

13 % worksheet for the specified row interval (s). Specify STARTROW... 

and 

14 % ENDROW as a pair of scalars or vectors of matching size for 

15 % dis-contiguous row intervals. To read to the end of the file ... 

specify an 

16 % ENDROW of inf.% 

17 % Example: 

is % [BoneLettering,InitialMajorAxisBoneplane,... 

InitialMinorAxisPerpendiculartoBone,FinalMajorAxis, ... 
FinalMinorAxis,DistancetoMark,NOTE] = importfile (' SWOTA BONE ... 
TABLE - OFFICE ( 1 ). xls ',' Sheet1 ', 2 , 256 ); 
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19 % 

20 % See also XLSREAD. 

21 

22 % Auto-generated by MATLAB on 2015/07/08 12:32:21 

23 

24 %% Input handling 

25 

26 % If no sheet is specified, read first sheet 

27 if nargin == 1 || isempty (sheetName) 

28 sheetName = 1; 

29 end 

30 

31 % If row start and end points are not specified, define defaults 

32 if nargin < 3 

33 startRow = 2; 

34 endRow = 25 6; 

35 end 

36 

37 %% Import the data 

38 [—i, raw] = xlsread (workbookFile, sheetName, sprintf ( ' A%d : G%d' , . . . 

startRow(1),endRow(1))); 

39 for block=2 : length (startRow) 

40 [—i, —i, tmpRawBlock] = xlsread (workbookFile, sheetName, ... 

sprintf ( 'A%d:G%d' ,startRow(block),endRow(block))); 

41 raw = [raw;tmpRawBlock]; %#ok<AGROW> 

42 end 

43 raw (cellfun (@ (x) — iisempty(x) && isnumeric (x) && isnan (x) , raw) ) = ... 

{”}; 

44 cellVectors = raw ( : , [1,7]); 

45 raw = raw ( : , [2,3,4,5,6]) ; 

46 

47 %% Create output variable 

48 I = cellfun (@ (x) ischar (x) , raw); 

49 raw (I) = {NaN}; 

so data = reshape ( [raw{ :}], size (raw) ) ; 

51 

52 %% Allocate imported array to column variable names 

53 BoneLettering = cellVectors (:, 1) ; 

54 InitialMa jorAxisBoneplane = data(:,l); 

55 InitialMinorAxisPerpendiculartoBone = data(:,2); 

56 FinalMa jorAxis = data(:,3); 

57 FinalMinorAxis = data(:,4); 

58 DistancetoMark = data(:,5); 

59 NOTE = cellVectors (:, 2 ) ; 


B.1.10 List loader 
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1 %% Import data from spreadsheet 

2 % Script for importing data from the following spreadsheet: 

3 % 

4 % Workbook: C:\Users\Ashkan\Google Drive\Code\Matlab\Bone ... 

drilling project\General\SWOTA BONE TABLE - OFFICE (1). xls 

5 % Worksheet: Sheetl 

6 % 

7 % To extend the code for use with different selected data or a ... 

different 

8 % spreadsheet, generate a function instead of a script. 

9 

io % Auto-generated by MATLAB on 2015/07/08 12:37:04 
n 

12 %% Import the data 

13 [-i, raw] = xlsread ( ' C : \Users\Ashkan\Google Drive\Code\Matlab\ . . . 

Bone drilling project\General\SWOTA BONE TABLE - OFFICE (1).... 
xls' ,' Sheetl' ); 

14 raw (cellfun (@ (x) -iisempty(x) && isnumeric (x) && isnan (x) , raw) ) = ... 

{”}; 

is cel IVe ct or s = raw (:, [1,7,8,9]); 

16 r aw = raw ( : , [2,3,4,5,6]) ; 

17 

is %% Replace non-numeric cells with NaN 

19 R = cellfun (@(x) —'isnumeric (x) && —dslogical (x) , raw) ; % Find non-... 

numeric cells 

20 raw (R) = {NaN}; % Replace non-numeric cells 

21 

22 %% Create output variable 

23 data = reshape ( [raw{ :}], size (raw) ) ; 

24 

25 %% Allocate imported array to column variable names 

26 BoneLettering = cellVectors ( : , 1) ; 

27 InitialMajorAxisBoneplane = data(:,l); 

28 InitialMinorAxisPerpendiculartoBone = data(:,2); 

29 FinalMa jorAxis = data(:,3); 

30 FinalMinorAxis = data(:,4); 

31 DistancetoMark = data(:,5); 

32 NOTE = cellVectors (:, 2 ) ; 

33 VarName8 = cellVectors (:, 3) ; 

34 VarName9 = cellVectors (:, 4) ; 

35 

36 %% Clear temporary variables 

37 clearvars data raw cellVectors R; 


B.1.11 Low pass filter function 2 
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1 function [y]=Low_filter_force(data) 

2 rng default; 

3 %x=data (:, 4 ); %noisy waveform 

4 x=data; 

5 h=fdesign.lowpass( 'Fp,Fst,Ap,Ast ',0.15,0.2,1,60); 

6 d=design(h, 'equiripple' ); %Lowpass FIR filter 

7 h=fdesign.lowpass('N,F3dB',12,0.0003); 

8 dl = design (h, 'butter ') ; 

9 y = filtfilt(dl.sosMatrix,dl.ScaleValues,x); 

10 %plot (x, ' b-. ' ) ; hold on; 

11 %plot (y, ' r ' , ' linewidth ' , 3) ; 


B.1.12 Labview oppener 


1 function [dataLV,H,MIN,SEC]=LTFO(filename) 

2 filelD = fopen (filename); 

3 fid=fileID; 

4 C = textscan (filelD, '%s %s'); 

5 Timeloc=ll; 

6 data_start=30; 

7 com=strfind(C{l,l}{Timeloc,l},' ,' ) ; 

8 colo=strfind(C{l,l}{Timeloc,1},': ' ) ; 

9 H= C{1,l}{Timeloc,1}((com+1):(colo(l,l)-l)); 

10 MIN= C{1,l}{Timeloc,1}( (colo (1,1)+1) : (colo (1,2)-1)); 

11 SEC= C{1,l}{Timeloc,1}((colo(1,2)+1):(colo(1,2)+8)); 

12 fid=fopen (filename) ; 

13 f id2=f open ( ' newf ile . txt ' , ' wt' ) ; 

14 id=0; 

is a=fgets (fid) ; 

16 while (ischar (a) ) 

17 id=id+l; 

is if id<24 

19 a=fgets (fid) ; 

20 continue 

21 else 

22 fprintf (fid2, a) ; 

23 end 

24 a=fgets (fid) ; 

25 end 

26 dataLV=csvread ( ' newfile . txt' ) ; 


B.1.13 Nan available correlation function 
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1 function r = nancorrcoef(x,y) 

2 %NANCORRCOEF - Compute Pearson's correlation for input vectors ... 

that contain NaN values 

3 % 

4 %Syntax: r = nancorrcoef(x,y); 

5 % 

6 %Inputs: X - First input vector 

7 % Y - Second input vector 

8 % 

9 %Output: R - Correlation coefficient 

10 % (Note: R is a scalar, not a 2 x 2 

u % matrix as in CORRCOEF) 

12 % 

13 %Example: x = randn(2 0,l); y = randn(2 0,l); 

14 % r = nancorrcoef (x, y) ; 

15 % 

16 %Other m-files required: none 

17 %Subfunctions: none 

is %mat-files required: none 

19 % 

20 %See also: CORRCOEF 

21 

22 %Author: Denis Gilbert, Ph.D., physical oceanography 

23 %Maurice Lamontagne Institute, Dept, of Fisheries and Oceans ... 

Canada 

24 %email: gilbertd@dfo-mpo.gc.ca Web: http://www.qc.dfo-mpo.gc.ca/... 

iml/ 

25 %September 2001; Last revision: 13-Sep-2001 by D. Gilbert 

26 

27 %Check number of input vectors 

28 if nargin y 2 

29 error ('Two input vectors are required'); 

30 end 

31 

32 %Conditions that must be satisfied to continue the calculations 

33 if size (x) y size (y) 

34 error (' Dimensions of input vectors must be the same !!!'); 

35 elseif (sum (isnan (x)) /length (x)) > 0.20 | (sum (isnan (y)) /length (y. .. 

)) > 0.20 

36 disp(' Input vectors have more than 20% NaN values...') 

37 disp( 'Verify if you could eliminate a few NaNs') 

38 disp('by interpolation or some other method') 

39 error ('Too many NaNs, processing stopped !!!'); 

40 elseif (sum (isnan (x)) /length (x)) > 0.05 | (sum (isnan (y)) /length (y. .. 

)) > 0.05 

41 warning (' Vectors have more than 5% NaN values, results may be... 

inaccurate !!!' ) 

42 end 

43 

44 %Only keep the common valid values in both input vectors 
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%Logical flag 


45 tf = -iisnan(x) & -lisnan (y) ; 

46 x = x (tf) ; y = y (tf) ; 

47 

48 meanx = mean(x); %Compute mean value of X 

49 meany = mean(y); %Compute mean value of Y 

50 

51 %Compute the arguments that go into the mathematical formula of R 

52 sx2 = sum ( (x-meanx) . "2) ; 

53 sy2 = sum ( (y-meany) . "2) ; 

54 sxy = sum ( (x-meanx) . * (y-meany) ) ; 

55 

56 % Mathematical definition of Pearson's product moment correlation... 

coefficient 

57 r = sxy . /sqrt (sx2*sy2) ; 


B.1.14 Omni data cutter 


1 clc 

2 clear all 

3 % [data_lv, h_lv,min.lv,sec.lv]=Dataanalyze('test.lvm'); 

4 %data.lv(: , 1)=data.lv(: , 1)+str2num(sec.lv)+str2num(min.lv)*60+.. . 

str2num(h.lv)*3600; 

5 data_omni=csvread ( 'R54.txt' ); 

6 points=[94169,119449,120426]; 

7 %data_omni(:,10)=12; 

8 data.omni(:,17)=data_omni(:,13)/1000+data.omni(:,12)+data_omni ... 

(: , 11)*60+data_omni(:,10)*3600; 

9 for q=l: length (data.omni) 

10 [x, y, z] =trans_omni ( [data.omni (q, 10) , data.omni (q, 11) , data.omni . . . 

(q,12)], [data.omni(q,1),data.omni(q, 2) , data.omni(q, 3)] , [ ... 
data.omni(q,7),data.omni(q,8),data.omni(q,9)],[-174 -81 ... 

0] ) ; 

u Location.drill(q, 1)=x;Location.drill(q, 2)=y;Location.drill(q. . . 

,3)=z; 

12 end 

13 

14 lenght_l_2= ( (Location.drill (points (1, 1) , 1) -Location.drill (points . . . 

(1.2) ,1))"2+(Location.drill(points(1,1),2)-Location.drill(... 
points(l,2),2))~2 +(Location.drill(points(1, 1) , 3)- . . . 
Location.drill(points(1,2) ,3))"2)".5 

15 lenght_l_2= ( (Location.drill (points (1,2) , 1) -Location.drill (points . . . 

(1.3) ,1))~2 +(Location.drill(points(1,2),2)-Location.drill(... 
points(l,3),2))~2 +(Location.drill(points(1,2) , 3)- . . . 
Location.drill(points(1,3) ,3)) ~ 2 )".5 

16 for num=points (1, 1) :points (1,2) 

17 v (num-points (1, 1) +1) = (data.omni (num, 4 ) " 2+data_omni (num, 5) "2+. . . 

data.omni(num, 6)~2) " . 5; 
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is end 

19 max_v=max(v) 

20 min_v=min (v) 


B.1.15 Plot generator function 

1 

%% force plot 





2 

d=[1 2 3 4 5 6 7 9 10 11 12 13 14 16 18 19 20 21 22 

23 

24 

25 

26 . . . 


27 28 29 30 31 32 34 35 36 37 38 39 40 41 43 44 

45 

46 

47 

49 ... 


50 51 52 53 54 55 56 57 60 61 62 63 65 66 67 68 

69 

70 

71 

73 ... 


75 76 77 78 79 80 82 83 84 85 86 87 90 91 92 93 

94 

95 

97 

98 ... 


99 101 102 103 104 105 016 107 108 110 111 113]; 





3 

for i=l: length (d) 





4 

x=Total_finaldata(d(i)).Time_lv_2ndcl iped; 





5 

y=Low_filter_force(Low_filter_force(Low_filter_ 

force ( 




Total_finaldata (d (i)) .Force_lv_2ndcliped))); 





6 

x=x-min (x); 





7 

y=y-min (y); 





8 

x=x. /max (x); 





9 

y=y. /max (y); 





10 

plot (x, y, 'b. ' ); 





11 

hold on 





12 

%title(num2str (d(i))) 





13 

%pause (1) 





14 

end 





15 

dd=[4 9 50 51 52 53 54 55 65 66 69 67 68 69 70 80 82 

83 

84 

85 

105] 

16 

for i=l:l%length (dd) 





17 

x=Total_f inaldata (dd (i) ) . Time_lv_2ndcliped; 





18 

y=Low_filter.force(Low_filter_force(Low_filter_ 

force ( 




Tot a l_f inaldata (dd (i) ) . Force_lv_2ndcliped) ) ) 

r 




19 

x=x-min(x); 





20 

y=y-min (y); 





21 

x=x. /max (x); 





22 

y=y./ max (y); 





23 

plot (x,y, 'r.' ); 





24 

hold on 





25 

%title(num2str(d(i))) 





26 

%pause(1) 





27 

end 






B.1.16 RPM function 1 
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1 

function [res]=rpm_final(w,t) 


2 

n=10000; 


3 

for i=l: length (w)/n 


4 

tim(i) =mean (t(((i—1)*n+l) 

((i)*n))) ; 

5 

ddd(i)=rpm2(w(((i—1)*n+l) 

((i)*n)))*60; 

6 

end 


7 

res=[ddd;tim]; 


B.1.17 RPM function 2 

1 

function [dd]=rpm2(data) 


2 

y=data-mean (data)/2; 


3 

%[y]=Low_filter_rpm(data) ; 


4 

d=0; 


5 

for i=l:( length (y)-1) 


6 

if ((y(i)*y(i+1))<0) 


7 

d=d+l; 


8 

end 


9 

end 


10 

dd=d/2; 


B.1.18 General SYNC function 


1 function [OUTPUT]=syncer(dataLVfA, haptic_A, clip_A) 

2 %OUTPUT=[Time, X,V, Force] ; 

3 %clip LV : 

4 newLV=dataLV_A(clip_A(1,1):clip_A(1,2),:); 

5 plot (dataLV_A(:,1),dataLV_A(:,2)); hold on ; plot (newLV(:,1),newLV. .. 

(:,2), '.r'); 

6 [—i, haptic_clip_l ] =min (sqrt (newLV (1,1) -haptic_A {:,!))); 

7 [-i, haptic_clip_2 ] =min (sqrt (newLV (length (newLV) , 1) -haptic_A (:,!))) . . . 

r 

8 new_data_LV=resample(newLV, haptic_clip_2-haptic_clip_l, clip_A. . . 

(1,2) - clip_A. (1, 1) ) ; 

9 %figure;plot(newLV (:, 1 ), newLV (:, 2 )); hold on; plot(new_data_LV ... 

(:,1), new.data_LV (:, 2) ,'. r '); 

10 size(new_data_LV); 

u size(haptic_A (haptic_clip_l:haptic_clip_2,2) ) ; 

12 size (haptic_A (haptic_clip_l : haptic_clip_2,3) ) ; 

13 OUTPUT= [haptic_A (haptic_clip_l : haptic_clip_2,7) , haptic_A ( . . . 

haptic_clip_l:haptic_clip_2,1) , haptic_A (haptic_clip_l: . . . 
haptic_clip_2,2),new_data_LV(:,2)]; 
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14 figure; 

15 plot(OUTPUT(:,2)*10); hold on; plot(OUTPUT(:,3)*10000,'r');plot(... 

OUTPUT(:,4), 'g' ) 


B.2 Augmented haptic project 
B.2.1 General code for type A 


1 clc 

2 clear all 

3 %close all 

4 %% 

5 [dataLV_A, H, MIN, SEC]=Labview2var( 'A_1.lvm' ) ; 

6 haptic_A=csvread('A_1.txt') ; 

7 %% 

8 [dataLV_A(:,2)]=Low_filter.force(dataLV_A(:,2)); 

9 dataLV_A (:,2) = (5601.7*dataLV_A (:,2)+6208.7)*0.00980665002864; 

10 m=f igure; plot (dataLV_A ( : , 1) , dataLV_A ( : , 2) ) ; 
u haptic_A (: , 3)=str2num(H); 

12 hapt ic_A ( : , 7 ) =hapt ic_A ( : , 3) *3600+hapt ic_A ( : , 4 ) * 60+haptic_A ( : , 5) + . . . 

haptic_A ( : , 6) *0.001; 

13 hold on; plot (hapt ic_A ( : , 7 ) , hapt ic_A ( : , 1) * 10, ' r ' ) ; plot (haptic_A . . . 

(:,7),haptic_A(:, 2)*10000, 'g') 

14 %% 

15 clip_A= [ 6.398e + 04 1.464e+05]; 

16 [m, OUTPUT] =syncer (dataLV_A, haptic_A, clip_A) ; 

17 saveas (m, 1 2 3 A_l_clip' , 'fig') 
is h=figure; 

19 s f_A_1 = fit([OUTPUT(:,2), OUTPUT(:,3)],OUTPUT(:,4),' poly23' ) 

20 plot (sf_A_l,[OUTPUT(:,2), OUTPUT (:, 3)], OUTPUT (:, 4 ) ) 

21 save ( ' A_1 ' , ' s f _A_1 ' ) 

22 saveas (h, ' A_1 ' , ' fig ' ) 

23 X=OUTPUT ( : , 2) ; 

24 dX=OUTPUT ( : , 3) ; 

25 F=OUTPUT ( : , 4) ; 


B.2.2 Force filter function 


1 d = designfilt('lowpassfir'); 

2 y = filtfilt(d,data); 

3 yl = filter (d,data); 
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4 subplot ( 2 , 1 , 1 ) 

5 plot ( [y yl]) 

6 title('Filtered Waveforms') 

7 legend (' Zero-phase FilteringConventional Filtering') 

8 

9 subplot ( 2 , 1 , 2 ) 

10 plot (data) 

u title('Original Waveform') 


B.2.3 Location filter function 


1 function [ data ] = fitrev( sf,x,y ) 

2 %UNTITLED5 Summary of this function goes here 

3 % Detailed explanation goes here 

4 data = sf.pOO + sf.pl0*x + sf.p01*y + sf.p20*x.~2 + sf.pll*x.*y +... 

sf.p02*y."2 + sf ,p21*x . "2 . *y 

5 + sf.pl 2 *x.*y ."2 + sf.p03*y.~3; 

6 end 


B.2.4 Labview opener 


1 function [dataLV,H,MIN,SEC]=Labview2var(filename) 

2 %The 6th culomn is the ABS time. The function creates a 'newfile. ... 

txt ' . 

3 filelD = fopen (filename); 

4 fid=fileID; 

5 C = textscan ( filelD, '%s %s' ); 

6 Timeloc=ll; 

7 data_start=30; 

8 com=strfind(C{l,l}{Timeloc,1},' ,' ); 

9 colo=strfind(C{l,l}{Timeloc,1},':'); 

10 H= C{ 1, l}{Timeloc, 1} ( (com+1) : (colo(l,l)-l)); 

11 MIN= C{1,l}{Timeloc,1}((colo(1,1)+1):(colo(1,2)-1)); 

12 SEC= C{1,l}{Timeloc,1}( (colo(1,2)+1) : (colo (1,2)+8)); 

13 fid=fopen (filename) ; 

14 fid2= fopen ( ' newfile . txt' , ' w ' ) ; 
is id=0; 

16 a=fgets (fid) ; 

17 while (ischar (a) ) 

is id=id+l; 

19 if id<24 

20 a=fgets (fid) ; 


88 









continue 


21 

22 else 

23 fprintf (fid 2 , a) ; 

24 end 

25 a=fgets (fid) ; 

26 end 

27 xx=csvread( 'newfile.txt' ); 

28 xy=xx(:, 1 )+( str 2 num (SEC) + 60 *str 2 num (MIN) + 3600 *str 2 num (H)); 

29 %Be sure the trend is accelerometer 0 ,accelerometer 1 , Voltage ... 

0 , Voltage 1 

30 %The structure is : ABS TIME ; Acceleration 0 ; Acceleration 1 ; ... 

Filtered 

31 %force ; RPM 

32 dataLV=[xy , xx ( :, 2 ) ] ; 

33 %dataLV=csvread ( ' newf ile . txt' ) ; 


B.2.5 File opener 


1 d = designfilt('lowpassfir'); 

2 y = filtfilt(d,data); 

3 yl = filter (d,data); 

4 subplot (2,1,1) 

5 plot ([y yl]) 

6 title('Filtered Waveforms') 

7 legend (' Zero-phase FilteringConventional Filtering') 

8 

9 subplot (2,1,2) 

10 plot (data) 

u title('Original Waveform') 


B.2.6 Calculator 


1 function [r_sq]=R_SQ_REG(input) 

2 %needs the data str made in clip_data_fix 

3 f orce=input .sf. pOO+input .sf.plO.* input .OUTPUT(:,2)+ input .sf.pO1.*... 

input .OUTPUT(:,3)+ input .sf.p2 0.*( input .OUTPUT(:,2) .~ 2 )+ input .... 
sf.pll.* input .OUTPUT(:,2).* input .OUTPUT(:,3) +input .sf.p02.* ... 
input .OUTPUT(:,3). ~2+input .sf.p21.* (input .OUTPUT( : , 2 ) . ~ 2 ) . * . . . 
input .OUTPUT(:,3)+ input .sf.pl2.* input .OUTPUT(:,2) .*( input .... 
OUTPUT(:,3) ."2) +input .sf.p03. *input .OUTPUT(: , 3) .~3; 

4 r_sq=( sum ( (force-mean ( input .OUTPUT(:,4))).~2) /sum (( input .OUTPUT ... 

(:,4) -mean ( input .OUTPUT ( :,4))) . " 2)); 
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B.2.7 Error calculator 


1 

A.sf.p00=mean 

( [A_l 

2 

B.sf .p00 =mean 

( [B_1 

3 

C.sf .p00 =mean 

( [C_l 

4 

D.sf .p00 =mean 

( [D_1 

5 

6 

A.sf ,p01=mean 

( [A_l 

7 

B.sf,p01 =mean 

( [B_1 

8 

C.sf,p01 =mean 

( [C_l 

9 

D.sf,p01 =mean 

( [ D_1 

10 



11 

A.sf ,pl0=mean 

( [A_l 

12 

B.sf ,pl0 =mean 

( [B_l 

13 

C.sf ,pl0 =mean 

( [C_l 

14 

D.sf ,pl0 =mean 

( [D_1 

15 



16 

A.sf ,p20=mean 

( [A_l 

17 

B.sf ,p20 =mean 

( [B_l 

18 

C.sf ,p20 =mean 

( [C-l 

19 

D.sf ,p20 =mean 

( [ D_1 

20 



21 



22 

A.sf,pll=mean 

( [A_l 

23 

B.sf,pll=mean 

( [B_l 

24 

C.sf,pll=mean 

( [C_l 

25 

D.sf ,pll=mean 

( [D_1 

26 



27 

A.sf ,p02=mean 

( [A_l 

28 

B.sf ,p02 =mean 

( [B_l 

29 

C.sf ,p02 =mean 

( [C-l 

30 

D.sf,p02 =mean 

( [D_l 

31 



32 

A.sf,p21=mean 

( [A_l 

33 

B.sf,p21 =mean 

( [B_l 

34 

C.sf,p21 =mean 

( [C-l 

35 

D.sf ,p21 =mean 

( [D-l 

36 



37 

A.sf ,pl2=mean 

( [A-l 

38 

B.sf ,pl2 =mean 

( [B_l 

39 

C.sf ,pl2 =mean 

( [C_l 

40 

D.sf ,pl2 =mean 

( [D-l 

41 



42 

A.sf .p03=mean 

( [A-l 

43 

B.sf.p03=mean 

( [B-l 

44 

C.sf.p03=mean 

( [C_l 

45 

D.sf.p03=mean 

( [D_l 

46 



47 




sf.pOO 

A-2 .sf.pOO 

A-3 

sf.pOO 

B_2 .sf.pOO 

B_3 

sf.pOO 

C_2 .sf.pOO 

C_3 

sf.pOO 

D_2.sf.pOO 

D_3 

sf.pOl 

A-2.sf.pOl 

A-3 

sf.pOl 

B_2 .sf.pOl 

B_3 

sf.pOl 

C_2 .sf.pOl 

C_3 

sf.pOl 

D_2 .sf.pOl 

D_3 

sf.plO 

A-2 . sf.plO 

A-3 

sf.plO 

B-2.sf.plO 

B-3 

sf.plO 

C_2.sf.plO 

C_3 

sf.plO 

D_2 . sf.plO 

D_3 

sf.p20 

A-2 . sf.p20 

A-3 

sf.p20 

B-2 . sf.p20 

B-3 

sf.p20 

C_2 . sf.p20 

C_3 

sf.p20 

D_2 . sf.p20 

D_3 


sf.pll 

A-2.sf.pll 

A-3 

sf.pll 

B_2.sf.pll 

B_3 

sf.pi1 

C_2.sf.pi1 

C_3 

sf.pi1 

D_2.sf.pi1 

D_3 

sf.p02 

A-2.sf.p02 

A-3 

sf.p02 

B-2.sf.p02 

B-3 

sf.p02 

C_2.sf.p02 

C_3 

sf.p02 

D_2.sf.p02 

D_3 

sf.p21 

A-2.sf.p21 

A-3 

sf.p21 

B-2.sf.p21 

B-3 

sf.p21 

C_2.sf.p21 

C_3 

sf.p21 

D_2.sf.p21 

D_3 

sf.pl2 

A_2.sf.pl2 

A_3 

sf.pl2 

B-2.sf.pl2 

B-3 

sf.pl2 

C_2.sf.pl2 

C_3 

sf.pl2 

D_2.sf.pl2 

D_3 

sf.p03 

A_2.sf.p03 

A-3 

sf.p03 

B-2.sf.p03 

B-3 

sf.p03 

C_2.sf.p03 

C_3 

sf.p03 

D_2.sf.p03 

D_3 


.sf.pOO A_4.sf.pOO]); 
.sf.pOO B_4.sf.pOO]); 
.sf.pOO C_4.sf.pOO]); 
.sf.pOO D_4.sf.pOO]); 

.sf.pOl A_4.sf.p01]); 
.sf.pOl B_4.sf.pOl]); 
.sf.pOl C_4.sf.pOl]); 
.sf.pOl D_4.sf.pOl]); 

.sf.p10 A_4.sf.plO]); 
.sf.plO B_4.sf.plO]); 
.sf.p10 C_4.sf.plO]); 
.sf.p10 D_4.sf.plO]); 

.sf.p20 A_4.sf.p20]); 
.sf.p20 B_4.sf.p20]); 
.sf.p20 C_4.sf.p20]); 
.sf.p20 D_4.sf.p20]); 


.sf.pll A_4.sf.pll]); 
.sf.pll B_4.sf.pll]); 
.sf.pll C_4.sf.pll]); 
.sf.pll D_4.sf.pll]); 

.sf.p02 A_4.sf.p02]); 
.sf.p02 B_4.sf.p02]); 
.sf.p02 C_4.sf.p02]); 
.sf.p02 D_4.sf.p02]); 

.sf,p21 A_4.sf.p21]); 
.sf,p21 B_4.sf.p21]); 
.sf.p21 C_4.sf.p21]) ; 
.sf.p21 D_4.sf.p21]); 

.sf.pl2 A_4.sf.pl2]); 
.sf.p12 B_4.sf.pl2]); 
.sf.p12 C_4.sf,pl2]); 
.sf,pl2 D_4.sf.pl2]); 

.sf.p03 A_4.sf.p03]); 
.sf.p03 B_4.sf.p03]); 
.sf.p03 C_4.sf.p03]); 
.sf.p03 D_4.sf.p03]); 
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48 A. OUTPUT 

49 B. OUTPUT 
so C. OUTPUT 
5i D. OUTPUT 


[ A_1.OUTPUT' 
[B_l.OUTPUT' 
[C_l.OUTPUT' 
[D_1.OUTPUT' 


A_2.OUTPUT' 
B_2.OUTPUT' 
C_2.OUTPUT' 
D_2.OUTPUT' 


A_3.OUTPUT' 
B_3.OUTPUT' 
C_3.OUTPUT' 
D_3.OUTPUT' 


_4 .OUTPUT' ] ' ; 
_4.OUTPUT']'; 
_4 .OUTPUT'] ' ; 
_4 .OUTPUT'] ' ; 
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APPENDIX C 


C++ CODES AND FUNCTIONS 


C.l Augmented Haptic 


1 /*The modified template code from open haptics library 

2 Library generated for haptic onmi device 2008 

3 */ 

4 #include <QHHeadersGLUT.h>/ /Include all necessary headers 

5 #include <HDU/hduMath.h> 

6 #ifdef _WIN64 

7 #pragma warning (disable:4996) 

8 #endif 

9 

io #include <cstdio> 
n #include <iostream> 

12 #include<f stream> 

13 using namespace std; 

14 

is #include <HL/hl.h> 

16 #include <HD/hd.h> 

17 

is #include <HDU/hduVector . h> 

19 #include <HDU/hduError . h> 

20 struct PointMass 

21 { 

22 

23 } ; 

24 hduVector3Dd m_position; 

25 hduVector3Dd m_velocity; 

26 HDdouble mjnass; 

27 HDdouble m_kStiffness; 

28 HDdouble imkDamping; 

29 hduVector3Dd springForce; 

30 hduVector3Dd damperForce; 

31 hduVector3Dd inertiaForce; 

32 hduVector3Dd v.location; 

33 

34 void HLCALLBACK buttonCB(HLenum event, HLuint object, HLenum ... 

thread. 
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35 HLcache *cache, void *userdata) ; 

36 / * .. . 

37 Servo loop thread callback. Computes a force effect. 

38 ****ic*1c*1c******ic*ic**jc*****1c-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-ki^-k-k-k-k-k-k-k-k, 

* / 

39 void HLCALLBACK computeForceCB(HDdouble force [ 3 ], HLcache *cache,... 

void *userdata) 

40 { 

41 


62 


63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 } 


PointMass *pPointMass = static_cast<PointMass *> (userdata); 

// Get the time a since the last update. 

HDdouble instRate; 

hdGetDoublev(HD_INSTANTANEOUS_UPDATE_RATE, SinstRate); 

HDdouble aT = 1.0 / instRate; 

// Set Variable Values 
m_mass = 0.01; 

hdGetDoublev(HD_NOMINAL_MAX_STIFFNESS, &m_kStiffness); 
m_kStiffness *= 0.5; 

m_kDamping = 2 * sqrt(m_mass * m_kStiffness); 

// Get the current proxy position from the state cache. 

// Note that the effect state cache is maintained in ... 
workspace coordinates, 

// so we don't need to do any transformations in using the ... 
proxy 

// position for computing forces. 
hduVector3Dd proxyPos; 

hlCacheGetDoublev(cache, HL_PROXY_POSITION, proxyPos); 

// Compute the inertia force based on pulling the point mass ... 
around 

// by a spring. 

//std::cout« m_kStif fness; 

springForce = m_kStiffness * (proxyPos - imposition); 

damperForce = -m_kDamping * mrvelocity; 
inertiaForce = springForce + damperForce; 

// Perform Euler integration of the point mass state. 
hduVector3Dd acceleration = inertiaForce / m_mass; 
mrvelocity += acceleration * aT; 
m_position += m_velocity * aT; 

// Send the opposing force to the device. 
force[0] += -inertiaForce[0]; 
force[1] += -inertiaForce[1]; 
force[2] += -inertiaForce[2]; 
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78 


79 

80 / * . . . 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


si Servo loop thread callback called when the effect is started. 

82 icicicifiricicifirif^icificific^icifiricicic-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-M-k-k-k-k-k-k-k-k. . . 

* / 

83 void HLCALLBACK startEffectCB(HLcache *cache, void *userdata) 

84 { 

85 fprintf(stdout, "Custom effect started\n"); 

86 

87 // Initialize the position of the mass to be at the proxy ... 

position. 

88 hlCacheGetDoublev (cache, HL_PROXY_POSITION, m_position) ; 

89 

90 m_velocity . set ( 0 , 0 , 0 ) ; 

91 } 

92 

93 

94 / * .. . 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


Servo loop thread callback called when the effect is stopped. 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
*/ 


kkkkkkkk , 


void HLCALLBACK stopEffectCB(HLcache *cache, void *userdata) 

{ 

fprintf(stdout, "Custom effect stopped\n"); 

} 


int main(int argc, char *argv[]) //initalizes the function "main" 

{ 

ofstream myfile; 
myfile.open ("example.txt"); 

HHD hHD; 

HHLRC hHLRC; 

//HDErrorlnfo error; 

hHD = hdlnitDevice(HD_DEFAULT_DEVICE); 

/*if (HD_DEVICE_ERROR(error = hdGetError())) 

{ 

hduPrintError(stderr, Serror, "Failed to initialize ... 
haptic device"); 

fprintf(stderr, "\nPress any key to quit.\n"); 
getch (); 
return -1; 
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120 


} / * * / 

hdMakeCurrentDevice(hHD); 


121 

122 

123 hHLRC = hlCreateContext (hHD) ; 

124 hlMakeCurrent (hHLRC) ; 

125 

126 HLuint effect = hlGenEffects (1) ; 

127 

128 

129 

130 /* Add a callback to handle button down in the collision ... 

thread. */ 

131 hlAddEvent Callback (HL_EVENT_1BUTT0ND0WN, HL_OB JECTLANY, ... 

HL_CLIENT_THREAD, 

132 buttonCB, 0 ); 

133 hlAddEventCallback (HL_EVENT_1BUTT0NUP , HL_OB JECT.ANY, ... 

HL_CLIENT_THREAD, 

134 buttonCB, 0); 

135 hlAddEventCallback (HL_EVENT_2BUTTONDOWN, HL_OB JECTLANY, ... 

HL_CLIENT_THREAD, 

136 buttonCB, 0); 

137 

138 hlBeginFrame () ; 

139 

140 hi Callback (HL_EFFECT_COMPUTE .FORCE, (HLcallbackProc) ... 

computeForceCB, 0); 

ui hlCallback (HL_EFFECT_START, (HLcallbackProc) startEffectCB, . 

0 ) ; 

142 hlCallback (HL_EFFECT_STOP, (HLcallbackProc) stopEffectCB, 0); 

143 

144 hlStartEf f ect (HL_EFFECT .CALLBACK, effect); 

145 

146 hlEndFrame ();/** / 

147 

148 fprintf (stdout, "Press any key to stop the effect\n"); 

149 //fprintf ("Press any key to stop the effect\n"); 

150 HDdouble key; 

151 key = 1 . 0 ; 

152 while (key = 1 ) 

153 { 

154 hduVector3Dd proxyPos; 

155 HLcache *cache; 

156 hlCacheGetDoublev (cache, HL_PROXY_POSITION, proxyPos); 

157 

158 std::cout « "proxy position = " « proxyPos[0] «", " « 

proxyPos [1] «", " « proxyPos [2] « " " « std::end. 

1 ; / * * / 

159 std::cout « "mass position = " « m_position [ 0 ] «", " . 

« m_position [ 1 ] «", " « m_position [2 ] « " " « ... 

std::endl; /**/ 
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160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 } 

181 

182 void HLCALLBACK buttonCB (HLenum event, HLuint object, HLenum ... 
thread. 


183 


HLcache *cache. 

void *userdata) 




184 { 







185 

if 

(event == HL_EVENT_1BUTT0ND0WN) 





186 

{ 






187 


/*fprintf(stdout, "Button 1 has 

been pressed.\n") 

r 



188 


float h; 





189 


h = 150.0; 





190 


//std::cout « "Angle 1 = " « pAngles->Thl « "\ 

n"; 



191 


//fprintf(stdout, "Angle 1 = %f\ 

n", pAngles->Thl) 

t 



192 


hduVector3Dd proxyPos; 





193 


hlCacheGetDoublev (cache, HLJPROXYJPOSITION, proxyPos) 

1 ; 


194 


std::cout « "Proxy x position = 

" « proxyPos[0] 

« 

1! 

\ n 



" « 





195 


"Proxy y position = 

" « proxyPos[1] 

« 

II 

\ n 



" « 





196 


"Proxy z position = 

" « proxyPos[2] 

« 

II 

II 



« std::endl;/**/ 




197 


//std::cout « "velocity = " « 

newVelFilt [ 0 ] «" 

II 

r 

« . 


newVelFilt [ 1 ] «", " « newVelFilt [ 2 ] « " " « std:.. 
:endl;/**/ 

198 //std::cout « "acceleration = " « accel[0] «", " « .. 

accel[l] «", " « accel[2] « " " « std: :endl;/**/ 

199 //std::cout « "force = " « inertiaForce [ 0 ] «", " « .. 

inertiaForce [ 1 ] «", " « inertiaForce [ 2 ] « " " « .. 
std::endl;/**/ 

200 //fprintf (stdout, "Stiffness = %f\n", m_kStif f ness) ; 


hlBegmFrame () ; 
hlCheckEvents(); 
hlEndFrame(); 

}; 

hlBeginFrame(); 
hlStopEffect(effect); 
hlEndFrame(); 

//fprintf(stdout, "Shutting down. . .\n"); 
//getch(); 

hlDeleteEffects (effect, 1); 

hlDeleteContext(hHLRC); 
hdDisableDevice(hHD); 

//return 0; 

//qhStartO; //Set everything in motion 
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201 

202 } 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 


219 

220 
221 
222 

223 

224 


225 

226 

227 

228 

229 

230 

231 


fprintf(stdout , "Mass = %f\n", m_mass); 

else if (event == HL_EVENT_1BUTT0NUP) 

{ 

} 

else if (event == HL_EVENT_2BUTT0ND0WN) 

{ 

/*fprintf(stdout, "Button 2 has been pressed.\n"); 
hduVector3Dd currentPos; 

hdGetDoublev(HD_CURRENT_POSITION, currentPos); 
std::cout « "Current x position = " « currentPos[0] « ... 
" mm\n" « 

"Current y position = " « currentPos[1] « ... 
" mm\n" « 

"Current z position = " « currentPos[2] « ... 
" mm" « std::endl; 

/*std::cout « "COM x position = " « COM[0] « " \n" « 

"COM y position = " « COM[1] « " \n" « 

"COM z position = " « COM[2] « " " « std:... 

:endl;/**/ 

//std::cout « "damping force = " « damperForce [ 0 ] «", ... 

" « damperForce [ 1 ] «", " « damperForce [2 ] « " " ... 
« std::endl;/**/ 

//std::cout « "spring force = " « springForce [ 0 ] «", "... 
« springForce [ 1 ] «", " « springForce [2 ] « " " «. . . 
std::endl;/**/ 
hduVector3Dd proxyPos; 

hlCacheGetDoublev(cache, HLJPROXY_POSITION, proxyPos); 

std::cout « "proxy position = " « proxyPos[0] «", " «... 

proxyPos [1] «", " « proxyPos [2] « " " « std:: end... 

1 ; / * * / 

std::cout « "mass position = " « imposition [ 0 ] «", " ... 

« imposition[1] «", " « imposition [2] « " " « ... 

std::endl; / **/ 

// basic file operations 
ofstream myfile; 
myfile.open ("example.txt"); 

myfile « "proxy position = " « proxyPos [0] «", " « . . . 

proxyPos [1] «", " « proxyPos [2] « " " « std mend... 

i; 

myfile.close (); 

} 
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New Mexico Tech 

SCIENCE • ENGINEERING - RESEARCH • UNIVERSITY 

Institutional Review Board 
for the Protection of Human Subjects in Research 

IRB Approval of Research Project 

REV. March 2013 PAGE 1 OF 2 

NMT's Institutional/Organization Number (assigned by DHHS): 

FWA00006606 

Valid Through: 12 July 2016 

IRB DATABASE No.: 

2014-10-004 


Title of research project: 

Robotic tools for bone drilling training and assessment 

Primary researcher name and title: 


David Grow 


Type of review sought by PI: 

^Initial Review of new research project ^Continuing Review for ongoing project 

Type of Initial Review performed: 

1X1 Full Board Review 

^Expedited Review - list IRB Members participating in review: 

Recommended Interval for Continuing 

Review: 

1 1 Every 6 months 

I^Every 12 months 

1 1 Further review not required unless 
protocol changes 

Research project qualifies for expedited IRB Review, according to the NMT IRB Policy, Appendix B, Category: 


^Collection of small samples of blood (<50 ml) 

I I Prospective collection of biological specimens for research purposes by non-invasive means 

^Collection of data through non-invasive procedures, including sensors, weighing, exercise testing, EEG, EKG, etc. 

^Research involving personally identifiable data, documents, or specimens collected solely for classroom/lab purposes 

^Collection of data from voice, video, digital, or image recordings made for research purposes 

1X1 Research on group characteristics or behaviors that also collects personally identifiable data 

I I Continuing review of research previously approved by the IRB: 

I I Where no additional subjects have been used and no additional risks have been identified 
Qwhere the remaining research activities are limited to data analysis 


IRB Review Findings: 

YES 

N/A 

1. The scientific questions addressed in this project have adequate merit to justify experimenting with human 
subjects. 


□ 

2. The risks to subjects are minimized and reasonable in relation to the anticipated benefits. 

13 

□ 

3. The selection of subjects is reasonable in relation to anticipated benefits. 

13 

□ 

4. The Pi's choice of consent documentation (e.g., written, oral, waived) is appropriate. 

13 

□ 

5. The project has appropriate protections to ensure the privacy of subjects and confidentiality of data. 

13 

□ 

6. No vulnerable populations are being used, OR there are appropriate safeguards for vulnerable populations. 

13 

□ 
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Phone: 575-835-5690 • FAX: 575-835-6934 
Email: jcarrillo@admin.nmt.edu • www.nmt.edu/human-subjects-research 























New Mexico Tech 

SCIENCE • ENGINEERING - RESEARCH • UNIVERSITY 

Institutional Review Board 
for the Protection of Human Subjects in Research 

IRB Approval of Research Project 

REV. March 2013 PAGE 2 OF 2 

NMT's Institutional/Organization Number (assigned by DHHS): 

FWA00006606 

Valid Through: 12 July 2016 

IRB DATABASE No.: 

2014-10-004 


Protocol has been reviewed by New Mexico Tech's IRB and is approved for 365 days from the date given below. 

Research must NOT continue after this date without continuing review and approval from the IRB. Research conducted without such review 
and approval runs the risk of being shut down and having all collected data destroyed. 

MONTH: October 

DAY: 24 

YEAR: 2014 

IRB Member (Printed! Jacob Carril'o 

Signature of IRB Member: 

Date Signed: 

IRB Member (Printed! ^rfrovlvn Zimrr _-rlv 

Signature of IRB Member: 

Date Signed: 
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Patient Code_ Date 

Surgical Skill Assessment Survey 


The Institutional Review Board for the Protection of Human Subjects at NM Tech has reviewed 
and approved this research project. By completing this survey, you agree to participate in this 
project. Your participation in this research project is strictly voluntary and you may choose not 
to participate by simply not completing the survey. You may also refuse to answer specific 
questions on this survey if you so desire. If you have any questions or concerns about this 
survey, please contact the researcher at 505-307-2918 or the New Mexico Tech IRB 
Administrator at 575-835-5690. 


Basic Information 

1. Gender: _ 

2. Age: _ 

3. Handedness: _ 

Previous Experience 

1. What is your current training level (e.g. 2 nd year resident)? 

2. Roughly, how many hours of training do you have in drilling real or synthetic bone? 


3. Do you have any prior experience working with haptic (force-feedback) robotic devices? If so, 


how many hours and what is this experience? 



ivttM rvirn TTm 












Subject Feedback 


1. Please rate the realism of the bone drilling experience compared to normally drilling in 


synthetic bone ( 1 [not at all realistic] - 10 [perfectly realistic] ): _(score) 


2. Please rate the realism of the bone drilling experience compared to normally drilling in human 


bone ( 1 [not at all realistic] - 10 [perfectly realistic] ): 


(score) 
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3. Please describe any sensations produced by the haptic (force-feedback) device that you found to 


significantly add to or detract from the realism of this training. 


Sample Recruiting Letter 

Dear Colleagues, 

We are currently conducting a human/machine interactive experiment in Dr. Grow's Human Robot 
Interfaces Laboratory. As part of this work, I will be conducting an experiment that investigates the use 
of robotic tools to measure skill in surgical drilling procedures and explore the use of these same tools 
for training of relevant skills. Results may lead to improved methods for assessing and training 
orthopaedic surgeons. 

This experiment will last 20-60 minutes. Only those who are 18 or over, and right-handed with no 
neurological disorders or injuries to their dominant hand are eligible to participate. If you are interested 
in being an experiment subject please contact me via email at pourkand@nmt.edu. 

Your participation in this experiment is entirely voluntary. If you have any questions or concerns, feel 
free to contact me or Dr. Grow at dgrow@nmt.edu or extension 5109. Your assistance in helping us 
meet our research goals would be greatly appreciated. 

Thank You, 

Ashlcan Pourkand 

IRB Protocol Number: XXXX 
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